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We show that D-dimensional de Sitter space is unstable to the nucleation of non-singular ge- 
ometries containing spacetime regions with different numbers of macroscopic dimensions, leading 
to a dynamical mechanism of compactification. These and other solutions to Einstein gravity with 
flux and a cosmological constant are constructed by performing a dimensional reduction under the 
assumption of g-dimensional spherical symmetry in the full D- dimensional geometry. In addition to 
the familiar black holes, black branes, and compactification solutions we identify a number of new 
geometries, some of which are completely non-singular. The dynamical compactification mechanism 
populates lower-dimensional vacua very differently from false vacuum eternal inflation, which oc- 
curs entirely within the context of four-dimensions. We outline the phenomenology of the nucleation 
rates, finding that the dimensionality of the vacuum plays a key role and that among vacua of the 
same dimensionality, the rate is highest for smaller values of the cosmological constant. We consider 
^L{, the cosmological constant problem and propose a novel model of slow-roll inflation that is triggered 

D ■ by the compactiflcation process. 

in 

■ I. INTRODUCTION 

Strong theoretical reasons for considering extra dimensions arise in many theories of physics beyond the Standard 
Model: the potential role of string theory as a theory of quantum gravity, the possible role of extra dimensions in 
addressing the hierarchy and flavor problems, and the ability of Kaluza-Klein theory to unify forces. Key to this 
O [ program is identifying ways in which extra dimensions can be hidden from our macroscopic view. The dominant 
^ ' paradigm of compactification has been greatly successful in producing reasonable four-dimensional effective theories, 
especially with recent developments in moduli stabilization in string theory. Indeed, there appears to be an embar- 
. rassment of riches, with compactifications using fluxes and branes leading to an enormous "landscape" of possible 
I> ' four-dimensional solutions 

, Determining the existence of phenomenologically acceptable four-dimensional vacuum states is important. Equally 
important is addressing the intrinsically cosmological question of the evolution into these and other vacua. How 
did our observable universe come to be effectively four-dimensional? If the laws of physics allow for more than one 
long-lived vacuum state, why do we find ourselves in a vacuum with a small cosmological constant and the particle 



\l ' spectrum of the Standard Model? These and other fundamental questions require an understanding of transitions 
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into different vacua, including states of different macroscopic dimensionality. 

In this paper we explore a new possibility for dynamical compactification: a spontaneous transition from a higher- 
dimensional de Sitter background to lower-dimensional spacetimes with vacuum energy that depends on details of 
the solution. The only ingredients in the simple models we consider are a positive cosmological constant in the 
original D-dimensional background space and a q-form gauge field strength. We show that this theory allows for the 
nucleation of regions containing a macroscopic D — g-dimensional spacetime with compactified g-spheres stabilized 
d by flux. This process is the inverse of the "spontaneous decompactification" studied by Giddings and Myers 0, and 
was first suggested therein. We explore the cosmology of these models, including a post-transition inflationary phase 
that could lead to a realistic universe. 

Other attempts to dynamically account for our observed four-dimensional world in a higher-dimensional theory 
include the suggestion of Brandenberger and Vafa [1] , who used the properties of intersecting strings to suggest that 
three spatial dimensions could naturally expand while others remained compact. Karch and Randall Q argued that 
a gas of branes with various dimensionalities would relax to one dominated by three-branes. Our proposal starts from 
a dramatically different initial state: empty de Sitter space in higher dimensions. Of course this is not necessarily a 
natural starting point from the perspective of string theory, which seems to favor Minkowski or anti-de Sitter vacua. 
Nevertheless, if some other process formed (for example) a six-dimensional de Sitter vacuum in string theory, our 
proposal could explain the transition from six to four dimensions and given the paucity of dynamical explanations for 
compactification, all such theories are potentially of interest. 

The post-transition solutions are derived from spherically symmetric black brane solutions of Einstein-Maxwell 
theory, and so share some properties with charged black hole solutions in four dimensions. In fact, we can explain 
most of the major elements of the _D-dimensional solutions and their nucleation by analogy with black holes. 

To begin, the interior of a black hole can be interpreted as a two-dimensional "big-crunch" cosmological spacetime. 
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The static metric of the Schwarzschild black hole inside of the horizon 



ds^ = - - 1^ dR^ + - 1 j dt' + R'dni (1) 

can be written in a suggestive form by making the change of coordinates 

X = t/AM, T = (IGM^ - MlRf^ , (2) 

and looking at the region just inside of the event horizon R < 2M (see also Ref. Q) 

ds^ ~ -dr^ + T^dx^ + AM'^dVil. (3) 

The r and x coordinates describe a two-dimensional cosmology, with the "extra dimensions" compactified on a two- 
sphere of momentarily constant radius. At t = the scale factor a = t goes to zero. In FRW cosmology this 
corresponds to the "big-bang," but in this case it is clear that the a = surface is merely a coordinate singularity 
coinciding with the location of the event horizon of the Schwarzschild black hole. The solution just outside the horizon 
can be obtained by analytically continuing r — > it, in which case r and x switch roles as the timelike or spacelike 
coordinates. In the region outside the horizon the radius of the two-spheres of course grows, "decompactifying" the 
extra dimensions as asymptotically flat four-dimensional space is reached far from the black hole. 
Exploring this scenario in more detail, we can write the metric Eq. [1] in the form 

ds^ ^ -dr^ + a^{T)dx^ + R'^{T)dnl, (4) 

which in the r, x plane describes a two-dimensional FRW universe in which the radius of the two-sphere R evolves 
with proper time r. The Einstein equations yield (see Ref. [6| for more details) 

i?+^ = -— = -^, a = R, (5) 
2R 2R dR ' ^ ' 

where it can be seen that the system of equations is that of a field R subject to the potential T4// = ^ logi? evolving 
in a two-dimensional FRW universe. Again, the "big-bang" at a = corresponds to an event horizon in the black 
hole geometry. The radius at which this occurs is in one-to-one correspondence with the mass of the black hole, 
and we can therefore classify solutions by stationary points in the motion of R (where, by the equations of motion, 
a — 0). As expected, the potential drives R ^ 0, corresponding to the singularity of the black hole. In addition to 
the non-singular big-bang surface, a non-singular big-crunch can be obtained by time reversing this solution. 

To generate the region of the black hole outside of the horizon, we use the intuition developed above, and analytically 
continue r — s- ir. This takes a ia and Veff — > —Veff in the equations of motion. Now, starting from a stationary 
point, R is pushed to infinity, corresponding to the asymptotically flat region far from the black hole. In this piecewise 
manner, it is possible to sew together the regions inside and outside the horizon by looking at the motion of R away 
from a stationary point in either the potential 14// or —Veff as shown in Fig,[TJ 

Within the context of black hole physics, we can find an example of the dynamical compactification mechanism 
mentioned above. In de Sitter space there is a non-zero probability to pair-produce charged black holes from the 
vacuum 0, [^, [l^, [ll|, [l^, [l^ • The two-dimensional FRW regions in the interior of each black hole are produced 
dynamically, and so black hole nucleation can be regarded as the dynamical compactification of two extra dimensions. 
The nucleation of many black holes with different charges leads to a fractally bifurcated structure of future infinity [3| 
reminiscent of eternal inflation (a process in which pocket universes of diverse properties are nucleated out of a 
background de Sitter space, dividing future infinity among many different vacua). However, since we are discussing 
black holes, the two-dimensional cosmologies always end in a big-crunch. 

A crucial new feature that arises when considering higher dimensional gravity is the existence of black brane 
solutions that are completely non-singular [15]. In contrast to the black hole solutions, this raises the possibility that 
there are true interpolating solutions between a region with D large dimensions and a lower dimensional non-singular 
cosmology. Semi-classical processes akin to black hole nucleation can then produce such regions dynamically. A 
sketch of the causal structure of the types of solutions we will discuss is shown in Fig. [21 The lower dimensional 
cosmology is contained within an event horizon, just as it was for the black hole solutions. Outside the horizon is an 
interpolating region to another event horizon that encloses an asymptotically _D-dimensional de Sitter region ^ . In this 



Because infinity is spacelike in de Sitter space, the approach to the asymptotically D-dimensional region is in this case in a timelike 
direction. The same will be true for the black hole geometries in de Sitter space, in the region outside of the cosmological horizon. 
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FIG. 1: A portion of the causal diagram for a Schwarzschild black hole. To the future of the event horizon, r is a timelike 
coordinate, and the r, x plane (which is shown here) resembles a big-crunch cosmology that begins from a completely non- 
singular "big-bang" at the event horizon. The value of R at the horizon is determined by the stationary point shown on the 
potential (filled circle). On subsequent surfaces of constant r, R evolves in the sketched potential V, eventually reaching a 
singularity as R is pushed to zero. To the past of the event horizon, r ir, and r becomes a spacelike coordinate. The 
value of R on surfaces of constant r is determined by its motion in the upside-down potential —V. This pushes J? oo as the 
asymptotically flat region far from the black hole is approached. 



D-q dimensions D dimensions 




FIG. 2: A sketch of the causal structure of a solution that interpolates between asymptotically D — q and 7?- dimensional regions 
across event horizons. On the left, a, D — q dimensional non-singular cosmological spacetime is located behind an event horizon. 
Outside this horizon is an interpolating region to another event horizon, that encloses an asymptotically _D- dimensional de Sitter 
region. The nucleation of such solutions from empty D-dimensional de Sitter space represents the dynamical compactification 
of some number of extra dimensions. 



geometry, either the D or p + 2-diniensional spaces can be accessed from the interpolating region. Similar proposals 
for embedding a lower dimensional cosmology within a higher dimensional space have been made in [l6l [TtI [TsI . [l^ . 

The lower-dimensional FRW universe can evolve to a variety of vacua in which the radion field is stabilized. Our 
model therefore describes a landscape that is simpler than string theory constructions, but nonetheless involves 
the existence of extra dimensions in a more fundamental way through the dynamical compactification mechanism. 
Using this mechanism, our model allows one to scan not only over different vacuum energies, but over fundamentally 
different spacetimes with different numbers of dimensions. This method of populating the states of a landscape is 
fundamentally different than the standard paradigm of four-dimensional eternal inflation. Finally, we note that the 
existence of landscapes should in fact be viewed as a generic consequence of theories with extra dimensions, underlying 
the importance of studying a number of different models. 

Summarizing the major elements of this picture, and their description in the remainder of the paper: 



4 



• It is possible to construct D-dimensional geometries piecewise by dimensionally reducing to a theory of lower- 
dimensional gravity with a radion field encoding the radius of some assumed number of spherically symmetric 
dimensions. The set of radion potentials define a landscape of lower-dimensional effective theories. The proper- 
ties of the dimensionally reduced theory are outlined in Sec. [TTl 

• A lower dimensional cosmological spacetime can be embedded behind an event horizon located in a higher 
dimensional geometry. The event horizons correspond to non-singular big-bang or big-crunch surfaces. We 
classify the possible ZJ-dimensional geometries under the assumption of an FRW metric ansatz in Sections llIIllIVl 
and|Vl We examine the case where the cosmology is homogenous but anisotropic in Appendix [X] Throughout 
we consider a zero or positive cosmological constant A > 0. A negative cosmological constant will yield solutions 
qualitatively similar to the A = case, and so we do not consider it here. 

• In a D-dimensional de Sitter space, there is a non-zero probability to nucleate solutions containing a non-singular 
lower-dimensional cosmological spacetime. We calculate the probabilities and describe the global structure of 
the D-dimensional spacetime after many nucleation events in Sec. I VII 

• In any theory that can access a number of sufficiently finely-spaced vacua, it is possible to address the cosmolog- 
ical constant problem. In Sec. lVIIi we briefly discuss the ability to obtain the observed value of the cosmological 
constant naturally from the landscape of radion potentials. 

• It is possible to embed an epoch of inflation inside of the lower-dimensional cosmological spacetime by adding a 
scalar field that obtains a negative mass squared inside an event horizon. This can be accomplished by coupling 
the scalar field to the I?-dimensional curvature and flux, as we describe in Sec. IVIIII 

II. DIMENSIONAL REDUCTION 

The theory that we consider throughout this paper is D-dimensional Einstein gravity (with and without a cosmo- 
logical constant) coupled to a g-form field strength. This system is described by the action 

S = ^ / d^.V^ (7t(^) - 2A - l^F,^) , (6) 

where Md is the D~dimensional planck mass (with SttGd = M'^^), and the tilde denotes that each quantity is 
evaluated in the D-dimensional Einstein frame. In each D-dimensional theory, there can be multiple g-form field 
strengths ranging from q — . . . , D — 2. 

Assuming g-dimensional spherical symmetry, we can write the metric in the product form 

ds^ = ~gP+^{K)dx^'dx'' + R^{x)dnl, (7) 

where 

D=p+2 + q, (8) 

dflg — d9\ + siv? OidOl -I- ... -I- sin^ 9i . . . siv? 9q-id9q, and ^ = 0, 1, . . .p -I- 1. We employ this decomposition into q 
and p + 2 because many of the objects we study are of codimension p (i.e. p-branes). The radius of the q-sphere i?(x) 
and p + 2-dimensional metric .g^ (x) will in general be a function of the p + 2-dimensional coordinates x. Magnetic 
g-form field strengths solving Maxwell's equations (trivially) and respecting the g-dimensional spherical symmetry are 
given by 

Fq = Q sin'^-Hi .. .sm9q^id9i ... A d9q. (9) 

For simplicity we do not consider electric fluxes in the following and construct solutions with a single fixed q (i.e. 
multiple g-form charges are not turned on simultaneously). 

Under our assumption of g-dimensional spherical symmetry, solutions to the full Z3-dimensional theory Eq. [S] can 
be described by the behavior of i?(x) in a p + 2-dimensional spacetime (in analogy with the discussion of black holes 
in the Introduction). In other words, we can perform a dimensional reduction by integrating over the angles of the 
g-sphcre. This yields a p + 2-dimensional theory of gravity coupled to a scalar field, the "radion," encoding the radius 
of the g-sphere at each point. This picture is crucial for describing how four-dimensional vacua and cosmologies can 
be obtained from a theory with extra dimensions. More generally, we will find that the dimensionally reduced theory 
is a powerful tool for constructing a variety of D-dimensional geometries not necessarily related to the description of 
four-dimensional universes. We now explore this picture in more detail, deriving the dimensionally reduced theory 
and describing its properties. 
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A. Defining the dimensionally reduced theory 



Because we have assumed (/-dimensional spherical symmetry, it is possible to recast the D-dimensional theory in 
the form of a p + 2-dimensional gravitational action coupled to the scalar field i?(x). We can write the _D-dimensional 
Ricci scalar and metric determinant as 



and 



(10) 



(11) 



The absence of a superscript on TZ and g indicates a p + 2-dimensional quantity. Integrating over the spherical 
coordinates and performing an integration by parts to remove the second derivative on R, the action Eq. [6] becomes 



s 



WR - q{q - l)R'^-^g''^{df,R)id^R) + q{q - _ 2Ai?9 



M 



2R1 



where ¥01(5') is the volume of a unit q-sphere 



Vol(S'' 



(12) 
(13) 



For p > 1, it is possible to perform a conformal transformation to the p + 2-dimensional Einstein frame 

5^1^ = {MoRfigt^u- 



(14) 



The absence of a tilde indicates that a quantity is evaluated in the p + 2-dimensional Einstein frame. Defining 

Mp+2 = Mc(Vol(5''))l/^ 



(15) 



the action becomes 



S 



+2, 



- '^'-^9^'^(d^R){d.R) - v{R) 



(16) 



where the potential for R in the Einstein frame is given by 



V{R) = ^ 



-q{q~l){MoRr^'-^ + ^^{MnRr^ + ^{MnR)-'^^^'n 



and we recover the standard Einstein-Hilbert action for the gravitational sector. 

It is possible to define a canonically normalized field (f> by making the change of variables 



(17) 



MdR — exp 



P 



q{p + q) Mp^ 



This yields the action 



S = 



2 



'Jl-^g'^-'{d,^)id,^)-Vi^) 



(18) 



(19) 



with 



F(0) 



-q{q~ l)exp ( -2 



p + q 



2A 

pq Mp+2 J ' 



exp ( —2 



p{p + q) Mp+2 



exp -2(p+ 1) 
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p{p + q) Mp+2 



(20) 



In the following sections, we will make frequent use of both the potential Eq. [2D]for and the potential Eq. [T7]for R. 
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I 1 \ \ / Increasing Q 



FIG. 3: The radion potential Eq. [20] with A = (left) and A > (right). In each plot, a number of potentials at fixed A are 
shown with successively larger values of Q from bottom to top. For A = the potential always has a negative minimum and 
approaches zero from below as ^ oo, while for A > there can exist a minimum of negative, zero, or positive energy and 
the potential approaches zero from above at large 0. 



The radion potential Eq. [^D] is defined in terms of a number of fixed parameters, such as the total number of 
dimensions D = q + p + 2 and the Z?-dimensional cosmological constant A. However, it also depends on the charge 
Q and the number of compact and non-compact dimensions q and p + 2 (under the constraint of fixed D), which are 
variable. This defines an indexed set of radion potentials, one for each allowed value of Q, p, and q. The potential 
Eq.[2Q]with fixed p and q is sketched in Fig.[3]for A = (left panel) and A > (right panel), in each case for a number 
of values of Q. Since it will be qualitatively similar to the A = case, we will not consider A < 0. The potential has 
the same qualitative features for all choices of p and q. 

Perhaps most interestingly, there are minima where the radius of the q-sphere can be stabilized, yielding p + 2- 
dimensional vacua. The set of vacua correspond to a "landscape" of lower dimensional theories ^. For A = 0, the 
vacua are always negative, decreasing in depth with increasing charge. For A > 0, the vacua are negative for small 
Q, but eventually reach zero and then become positive as Q increases. For large enough Q, the minimum disappears 
completely. The existence of 4-dimensional vacua with positive vacuum energy will be important in our discussion of 
cosmological solutions. 

The locations of potential extrema can be determined most easily from Eq. [T71 We will refer to the location of 
potential minima as i?_ or 0_ and potential maxima as or (j)^ . Setting the derivative of the potential with respect 
to R equal to zero, the extrema are given by the roots of 



B. A landscape from radion potentials 



A = ^(P + q){q - l)R-' ~\{P+ l)Q2M^2('-l)i?-2«. 



(21) 



In Fig. m the right hand side of this equality is plotted for various Q. 

When A = 0, there is a single extremum, whose position can be found directly from Eq. [^T] 



This is a global minimum of the potential. Inserting this back into Eq. 1171 the value of the potential is 




(22) 




(■i+p) 



(23) 



^ The existence of a landscape of vacua is strongly associated with string theory. However, a landscape of vacua can be viewed as a generic 
consequence of theories that have extra dimensions. There is even a landscape of two-dimensional vacua in the completely pedestrian 
four-dimensional Einstein-Maxwell system. 
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FIG. 4: Extrema of the radion potential are located at the intersection of each curve (Eq. [2T]for {p — 2, q — 2}) with a line of 
fixed A. There is only one intersection at A = and zero, one, or two intersections at A > 0. When there are two intersections, 
the first corresponds to a minimum (referred to as R- or 0_) and the second to a maximum (referred to as R+ or Note 
that the location of the first intersection is relatively independent of A. 



The depth of the potential minimum decreases with increasing Q, as can also be seen in Fig. |31 

For fixed A > 0, there are zero, one, or two roots depending on the value of Q. This is seen in Fig d] by locating 
the intersections of the curve with lines of fixed A. The maximum value of Q for which a root exists at fixed A can 
be found by determining the maximum of Eq. 1211 There are two roots as long as 

2l-^/^(g-l)^-V^(p + g)^/^ /Af^y^^^)/^ 



g«/2(H-p)i/2 \^ A 

which grows for decreasing A and increasing p and q. 

When it exists, the location of the smaller root is not strongly affected by the value of A, and given approximately by 
Eg. 1221 This is a local minimum of the potential. The location of the second root, when it exists, can be approximated 
by neglecting the Q-dependent term in Eq. [51] (since the second root is at relatively large R, this term is subdominant), 
yielding 



{p + q)iq-l) 



Inserting this into Eq. 1171 the height of the potential maximum is given by 



In addition, this roughly bounds the possible vacuum energy. As we will discuss in Sec. IVI) this bound ensures that 
the entropy associated with solutions at the positive potential minima is always greater than the entropy associated 
with a D-dimensional de Sitter space with cosmological constant A. 

Another property of note is the asymptotic behavior of the potentials as — > oo. For A = the negative term in 
the potential dominates at large cf), and the potential approaches zero from below. In the solutions of the following 
sections, we will see that this implies that the radius of the g-spheres increase in a spacelike direction, consistent with 
the assumption of asymptotically flat space. 

For A > 0, the term in the potential proportional to A dominates at large (/), and the potential approaches zero 
from above. We will find that this implies the radius of g-spheres increase in a timelike direction, consistent with a 
ZJ-dimensional de Sitter space. 

Therefore, both the vacuum structure and asymptotic properties of the potentials will play prominent roles in the 
construction of solutions in the dimensionally reduced theory, to which we now turn. 
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III. SOLUTIONS WITH A p + 2-DIMENSIONAL FRW METRIC ANSATZ 



We generate solutions by assuming that the p + 2-dimensional spacetime can be foUated into either spacehke 
or timehke p + l-dimensional surfaces of homogeneity. A fuU classification of homogenous solutions in arbitrary 
dimensions is beyond the scope of the present paper and we will restrict ourselves to two different classes of metrics. 
In the following three sections, we will treat the FRW case, classifying all possible solutions for open, closed, and flat 
spatial sections. In Appendix \^ we treat the case of an anisotropic metric, assuming that the spatial sections are 
flat (Bianchi Type I in four dimensions). The set of solutions generated from these two different forms of the metric 
will include both known extremal and non-extremal charged p-brane solutions, as well as a number of new solutions 
for A = and A > (again, the solutions for A < will be qualitatively similar to those for A = and we do not 
consider them in detail here). 

The evolution of the radion field on successive p + l-dimensional surfaces specifies the location and properties of 
singularities, event horizons, and asymptotic infinity. Specializing to the FRW metric ansatz, in addition to the radion 
field, there is one scale factor describing the geometry. We define the relevant equations of motion in Sec. IIII Al 

As illustrated in the introduction, regions with spacelike and timelike surfaces of homogeneity can be matched across 
event horizons in the full geometry. From the perspective of the dimensionally reduced theory, these are big-bang and 
big-crunch surfaces of p + 2-dimensional cosmological spacetimes. We describe the general circumstances in which 
this correspondence is valid in Sec. IIII Bl 

With the relevant formalism in hand, we will move on to discuss the particulars of the solutions in Sec. lIVI andlVl 
Although it is in most cases impossible to find a closed analytic form for the metric, we numerically generate examples 
of various classes of solutions. See Ref. "201 for a similar numerical study of the evolution of the radion in the presence 
of an FRW metric ansatz. The solutions relevant for the dynamical compactification mechanism are discussed in 
Sec. |Vl although the more simple geometries of Sec. I IVI provide a nice illustration of our method of finding solutions. 

In Sec. lIVI we catalog solutions with A = 0. These include the x 5"' compactification solutions (corresponding 

to sitting at minima of the radion potential) as well as various extremal and non-extremal p-branc solutions in 
asymptotically flat ZJ-dimensional space. We illustrate in detail the construction of the p-brane solutions from the 
dimensionally reduced theory. In addition to the known solutions, we find a novel non-singular geometry in which 
two identical asymptotically flat regions are connected across a pair of event horizons. 

Sec. |V] treats the case of a positive Z?-dimensional cosmological constant A. The compactification solutions in this 
case can be either AdSp+2 x S"*, M.p+2 x S"*, or ^5*^+2 x S"', corresponding to the negative, zero, or positive vacuum 
energy minima and the positive energy maximum of the radion potential. Among the geometries that contain an 
asymptotically D-dimensional de Sitter region will be the "interpolating" solutions mentioned in the introduction. 
These completely non-singular geometries interpolate across event horizons between an asymptotically D-dimcnsional 
de Sitter space and a,p-\- 2-dimensional open FRW universe that can evolve to dSp+2 x S'^. This geometry will form 
the basis of the dynamical compactification mechanism described in Sec. IVII 



A. Equations of motion 

For spacelike surfaces of homogeneity, the metric is that of an FRW universe 

ds^ = -dr^ + a{Tf [dx^ + Sl{x)dnl] , (27) 

where 

Si = {sinx, X, sinhx}, (28) 

for closed, flat, or open spatial slices respectively. The resulting equations of motion for the field 0(r) and the scale 
factor a(T) are given by 



^+(p+l)«0=_Af2;fV-', (29) 



where k — {—1,0,+!} for open, flat, and closed slicing respectively and V{(t)) is given by Eq. [20l We will also have 
occasion to use the second order equation for a, given by 
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We refer to surfaces where a = as big-bang or big-crunch surfaces depending on their orientation in time. 

When the surfaces of homogeneity are timehke, r becomes a spaceUke variable. The metric in these case can 
be obtained by analytic continuation. Closed FRW universes always yield spaces with closed timelike curves upon 
analytic continuation, and we do not consider them here. Starting with an open FRW universe, and taking 

yields 

ds^ = -a{T)^dx'^ + dr^ + a^T)"^ cosh^ xd^l, (33) 

where the surfaces of homogeneity are hyperboloids with spacelike norm. The equations of motion for the scale factor 
a(r) and (/)(t) are given by Eqns. [29l [30l and [31] with V —V. Recall that a similar flip of the potential occurred 
in the black hole example. We will refer to this as the "upside-down" or "Euclidean" potential. 

In the case of a flat universe, the analytic continuation t ir also takes V ^ —V in the equations of motion. 
Continuing one of the cartesian coordinates z — > iz as well, the metric becomes 

ds^ = -a{T fdx^ + dr^ + a(r)\2 {^91 + sinh^ OidVLp^i) . (34) 

The surfaces of homogeneity in this case are planes with spacelike norm. 

Solutions can be generated by numerically evolving the equations of motion. To do so, it is useful to define the 
dimensionless variables 

A, . V _ ^(P+2)/2 _ ^(P+2)/2 



Mlil ' Mli\ 



in terms of which the equations of motion become 



a 



^+{p+l)^^^~V', (36) 
a 



4' (37) 




a 



a pip + 1) 



-p^^ + 2V]. (38) 



B. Non-singular big-bang and big-crunch surfaces 

Just as for the black hole, we can match different regions across non-singular big-bang or big-crunch surfaces, which 
from the perspective of the full D-dimensional spacetime correspond to event horizons. However, unlike the example 
of the black hole, the equations of motion no longer guarantee that stationary points in the radion field occur when 
a = 0, and we must worry about the existence of curvature singularities. We will again find that a non-singular 
initial value surface is only possible by requiring a stationary point (f) — at a = 0. In addition, the scale factor 
must have the universal behavior a = t in the case of an open FRW cosmology or a = exp(i/T) in the case of a flat 
FRW cosmology. In each of these cases, the initial value surface is null and corresponds to an event horizon when the 
solutions are lifted to £)-dimensions. 

To determine the location of a singularity, we can write the Ricci scalar in terms of the fleld (f) as 

Therefore, as long as the kinetic and potential energy of the field remain bounded, there is no curvature singularity. 

A trivial way to satisfy this condition is for the field to sit at a critical point of the potential for all times, yielding 
a p + 2-dimensional de Sitter, Minkowski, or Anti-de Sitter space depending on the sign of the potential in regions 
where r is timelike. Such solutions were first discussed in Ref. 21], and are sometimes referred to as Freund- Rubin 
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compactifications. These spaces are homogenous in time as well as along the spatial sections, and so admit multiple 
time slicings, corresponding to the choices for spatial curvature in the FRW metric ansatz. de Sitter space admits a 
closed, flat, or open foliation while Minkowsi and Anti-de Sitter admit a flat or open foliation only. In each of these 
cases, the a = surface is simply a coordinate singularity which can be continued across. 

There are also non-singular solutions when the field has non-trivial dynamics. It is necessary to examine each 
ansatz for spatial curvature separately when evaluating these solutions, since the curvature term in Eg . 1301 tvpicallv 
dominates as a — > 0. 

For a closed universe (k = +1), there are solutions where a ^ only if the kinetic or potential energy of the 
field diverge faster than +l/a^ (otherwise d/o becomes imaginary). Since the energy density in the field becomes 
unbounded, such solutions produce a curvature singularity at a = 0. Otherwise, the scale factor has a non-zero 
minimal value, going from contracting to expanding and producing a "bouncing" cosmology with no a = surfaces. 
We will not consider this case further. 

In an open universe, the kinetic and potential energy of the field can remain finite as a — > 0. In fact, it is possible 
to have curvature alone, in which case the Milne patch of Minkowski space is generated. Here, the a = surface 
is obviously non-singular, and simply corresponds to a boundary of the coordinate patch. Putting (p back in, the 
absence of a curvature singularity requires an everywhere-finite energy density in the field. This in turn guarantees 
that the curvature term in Eq. [30] will always dominate as a — s- 0, yielding a universal behavior for the scale factor 
a = T as T —^ 0. Inserting this into the equation of motion Eq. [21] for 0, consistency requires that (f){T — 0) — due 
the divergent friction/anti-friction coefficient at a = 0. Solutions with these boundary conditions have a completely 
non-singular a = surface. 

In a flat universe, it is again necessary to set (j) = when a = to avoid a divergent friction/anti- friction coefficient. 
From Eq. 1301 this requires that y > to avoid an imaginary Hubble parameter. These conditions imply a universal 
behavior near the big bang surface since the positive vacuum energy comes to dominate at early enough times. In 
this case, one obtains the flat slicing metric coefficient a = exp^Hr) as r — > — oo. This can be recognized as the 
metric coefficient for de Sitter space in the flat slicing. The t — oo surface is merely a boundary of the coordinates, 
but in contrast to the case of negative spatial curvature, the initial conditions surface is pushed infinitely far to the 
past T — > — cx) along the worldlines of comoving observers (although boosted observers will encounter the initial value 
surface in finite proper time). 

This is an exhaustive list of the possible non-singular initial value surfaces under the FRW assumption for the 
p + 2-dimensional metric. 

IV. SOLUTIONS FOR A = 

In this section, we focus on geometries in the absence of a D-dimensional cosmological constant. In addition to 
the AdSpj^2 X S"^ compactification solutions, there are three types of geometries that include an asymptotically flat 
Z?-dimensional region: singular solutions (timelike and spacelike), oscillatory solutions, and extremal solutions. We 
will construct examples of each of these geometries in a piecewise manner by matching solutions for a(r) and 4>{t) 
across non-singular big-bang and big-crunch surfaces (event horizons) where = and a — 0. A representative 
example of the trajectory </)(t) in the presence of the potential Eq. [5n]is sketched in Fig. [5] for each type of solution. 
The causal structure of the corresponding ZJ-dimensional spacetimes are shown in Fig. [5] Each of these solutions is 
an extended object of codimension p (i.e. a p-brane) embedded in D-dimensional flat space. 

A few properties of these solutions can be understood without a detailed analysis. First, as described in Sec. IIIBi 
when A = 0, V{(f)) always has a negative global minimum and approaches zero from below as (f) oo. Recall 
that the motion is in the potential V{(j)) in regions where r is timelike and — T^(0) in regions where r is spacelike. 
Therefore, the radius of the g-sphere will always approach Z3-dimensional space in a spacelike manner, consistent with 
an asymptotically flat solution. If the field energy diverges along the trajectory of (p, the orientation of the potential 
also determines if the resulting curvature singularity is spacelike or timelike. In geometries without an event horizon 
(where there is no stationary point in </>), the naked singularities are timelike, while singularities cloaked behind an 
event horizon are necessarily spacelike. 

We begin our detailed discussion of the various geometries with the compactification solution, where the field sits 
at the potential minimum for all r. Since the minimum is negative, the full solution is a p -I- 2-dimensional Anti-de 
Sitter space with a stabilized g-sphere of constant radius: AdSp+2 x <5"'- The scale factor can be found from Eq. [311 
which if T is timelike, becomes 



(40) 
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FIG. 5: A sketch of the solutions for A = with a flat or open FRW ansatz for the p + 2-dimensional metric. Beginning 
from stationary points 4> = (filled circles), the evolution is in the potential V for r timelike (trajectories above the potential) 
or —V for r spacelike (trajectories below the potential). Singularities along a trajectory are denoted by an asterisk. The 
compactification solution sits at the potential minimum for all r. The singular (timelike and spacelike), extremal, and oscillatory 
solutions all have a region where r is spacelike and oo, universally approaching the attractor solution Eq. 1531 Continuing 
across an arbitrary stationary point to a region where r is timelike will result in a spacelike singularity. There are special 
stationary points from which oscillates about the potential minimum (indicated by the double arrows on the oscillatory 
trajectory). For a flat FRW metric ansatz, there are solutions where r is everywhere spacelike, and the field interpolates 
between the potential minimum and (f> ±cx). Solutions without a stationary point are always singular (timelike). 



oscillatory 




extremal (p odd) extremal (p even) 



FIG. 6: The D-dimensional causal structure of the solutions discussed in the text and shown in Fig. [5] In the top row, we 
show the timelike singular (left), spacelike singular (center), and oscillatory (right) solutions generated from an open FRW 
metric ansatz. In the bottom row are the extremal solutions generated from a flat FRW metric ansatz for p odd (left) and p 
even (right). Event horizons correspond to the non-singular big-bang and big-crunch surfaces at the stationary points of the 
p + 2 dimensional cosmological solutions, and are color-coded to correspond to the stationary points in Fig. [S] Singularities 
occur when the energy density in (j> diverges at a = 0, while the asymptotically flat _D- dimensional regions are described by the 
attractor solutions at large (f). 
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where 

For an open FRW metric ansatz, the solution is 

a(T) = Lj-^ sin (cjaT) . (42) 

There is no sohition for a flat or closed FRW ansatz with r timelike (the Hubble constant would be imaginary). 
However, for spacelike r there are solutions, corresponding to different time slicings of the AdS space. 

In all other cases, the field has dynamics. As discussed in Sec. IHI Bl the behavior of the field will be quite different 
for a closed, fiat, or open FRW metric ansatz. We will only consider the open and fiat cases in detail since they allow 
for the possibility of non-singular big-bang and/or big-crunch surfaces. 



A. Open FRW metric ansatz 



We first specialize to the case of open surfaces of homogeneity. At large </>, far from the potential minimum, consider 
evolution from a generic set of initial conditions {a(Ti) , (f>{Ti) , (j){Ti)} . To obtain asymptotically flat D-dimensional 
space, r must be a spacelike variable, and the evolution is in the inverted potential. Evolving to either earlier or later 
T from Ti, the field will overshoot the inverted-potential maximum if is to large. This causes the field to run off to 
(/) — > ~oo, resulting in a curvature singularity. An example of such a "singular (timelike)" solution is shown in Fig. [51 
and the resulting D-dimensional spacetime is shown in the top left panel of Fig. [6] (which, as could be ascertained 
from the naming scheme, has a naked timelike singularity). 

However, as described in Sec. IHI Bl these singularities can be avoided by ensuring that (j) — when a ~ 0. Specifying 
boundary conditions for cf) and a at the position of a stationary point in the field ((/> = 0), we can pick out the non- 
singular solutions by setting a — 0. If the stationary point occurs at some </) > then the field will evolve to 4> oo, 
and therefore all such solutions have an asymptotically Z?-dimensional region containing an event horizon. At the 
event horizon, it is possible to continue to a region where r is timelike, and the motion is in the original potential. 
By analyzing the motion away from stationary points at various cf) for both orientations of the potential, it is possible 
to build up the full D-dimcnsional geometry piece by piece. A similar classification of solutions by their stationary 
points was made in the black hole example described in the Introduction. 

All of the solutions will have a region with spacelike r that evolves towards </> ^ oo. At large (j), the radion potential 
Eq. [201 becomes dominated by a single term 

y^,,._, «'M-^K .,UMIj^). (43) 

2 \ \ pq Mp+2j 

It is well known that there are attractor solutions for the evolution of a scalar field in the presence of an exponential 
potential [2^ . [23| . Generic initial conditions will approach these fixed points in phase space, and the asymptotics can 
be determined. Each geometry should therefore have the same behavior in the asymptotically fiat region. 

We will find it convenient to work in terms of the dimensionless set of FRW variables defined in Eq. [35l and begin 
by making the ansatz that the scale factor has power law behavior (see e.g. 0, [13]) 

a = ao(cf )*. (44) 

For an exponential potential of the form we have been studying 

V{$) = Ae-'''!', (45) 

the field equation is satisfied for 

0=-log(cf), (46) 
7 



with the constant c fixed to be 
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Since there is both curvature and field energy, it is necessary to determine which will dominate the dynamics of the 
scale factor. Substituting for (f) and a from above, we can write Eg. 1301 as 



2 



4(5 1 



(48) 



where the first term represents the contribution from the field and the second from the curvature term. It can be seen 
that the field dominates the dynamics of the scale factor ior 6 > 1. If this is the case, using Eq. 231 S is found to be 

If this is less than one for a given potential, then neglecting the curvature term is not self consistent, and the equations 
of motion imply that S — 1. This analysis was performed by Ref. [S], who also showed numerically that these two 
behaviors are attractors in phase space. 

Setting — Mp_|_2M|,/2 and comparing with the analysis above, we have 



A^q{q~l), j^2jP±^, =^<1, (50) 

pq (p7^) p + q 



implying that the attractor solution has 6 = 1. The constant c is given by 



c= f'^f-'\ (51) 

From Eq. [181 we can find R{t) 

MoRir) = (cf)^^ . (52) 
Using this, we can find the full metric Eq. [Tlin the D-dimensional Einstein frame 

M^ds^ = (cf)-2?T? [df^ + a2^2-2 ^_^^2 _^ ^Qgj^2 ^dnl)] + {cff^dn^. (53) 
Calculating the D-dimensional Ricci scalar, in all cases it is proportional to 

7?,(^) cx f~Ffe cx 4^, (54) 

and the metric is asymptotically flat as f — > cx) as expected. 

Returning to the stationary point = 0, a = of the non-singular solutions (i.e. the event horizon in the D- 
dimensional geometry), we can analytically continue to timelike r. Here, there is no universal behavior, and different 
stationary points will produce qualitatively different solutions. For most choices, the solutions become singular. 
These portions of the "singular (spacelike)" solutions describe a region containing a spacelike singularity separated 
from asymptotically flat Z?-dimensional space by an event horizon. For some special choices of the stationary point, 
completely non-singular solutions can be generated. These "oscillatory" solutions describe a region between two event 
horizons in which the radius of the q-sphere oscillates. Across each event horizon is a different region of asymptotically 
flat D-dimensional space. We now analyze these cases in more detail. 

For stationary points in the near vicinity of the negative potential minimum, Eq. [5T] implies that the scale factor 
a is bounded, and possess two zeroes. For small amplitude motion of (p, then \V\ ^ cfP' , and the solution is given 
approximately by Eq.|42l Substituting Eq. [42] into the field equation Eg . [29] and Taylor expanding the potential about 
the minimum we obtain 

ij>+{p+l) COi{L0aT)^ = - Afl^f |F"(0_ )U-<P-). (55) 

This has an exact solution in terms of Gegenbauer polynomials 

cj,~(P^= ACl/^ [cos {uJaT)], (56) 

where A is the amplitude of oscillations and the index a is determined by the positive root of 



(57) 
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FIG. 7: The oscillatory solutions for timelike t with A = 0, Q = 1, p = 5, and q = 2. The solutions for a(f ) (left in each row) 
and (j>{f) (right) describe the spacetime between two event horizons of the D-dimensional geometry (see Fig.[S|. Event horizons 
are located at the zeros in a. The top panel shows the oscillatory solutions with 4 nodes, generated from the initial condition 
00 — .105. The bottom panel shows the oscillatory solution with 2 nodes generated from the initial condition 4>o ~ .523. A 
scan over initial conditions reveals that these are the unique non-singular solutions. 



The ratio of the potential to its second derivative can be evaluated using the results of Sec. Ill Bl yielding 

\V"{^^)\ 4(p+l) 



(58) 



\V{<P-)\ P ' 

Finally, solving for cr, we find 

a = p+l. (59) 

Amazingly, for the radion potentials a is always an integer. For a integer, the Gegenbauer polynomial has a nodes 
and is zero at the endpoints of the motion in a. Such solutions are completely non-singular. 

Considering finite size excursions away from the minimum, the corrections to the frequency of oscillations of (f> 
become important. As a approaches its second zero, (j) is no longer zero, and the (p term in the field equation Eq. [29] 
yields a divergent anti-friction. This corresponds to the appearance of a spacelike curvature singularity (since we are 
evolving in a region where r is timelike) in the geometry. A non-singular stationary point to the right of the potential 
minimum therefore generically connects the attractor solutions discussed above to a region containing a spacelike 
singularity. The Z)-dimensional causal structure of such singular (spacelike) solutions is shown in the top center panel 
of Fig. H 

Although a generic choice of initial conditions leads to a singularity, there are special choices that lead to a completely 
non-singular evolution of the field. For an initial condition to the right of the minimum, the period of oscillation in 
(j) increases monotonically with increasing amplitude (the potential becomes more gradual at large 0). By increasing 
the period of oscillation, eventually another harmonic of the period in a will be reached. The asymmetry of the 
potential forbids solutions with an odd number of nodes, and it follows that there will be initial conditions leading to 
non-singular solutions with a = 2, 4, . . .p — 1 for p odd and cr = 2, 4, . . .p for p even. These non-singular "oscillatory" 
solutions connect two asymptotic regions where (j) oo across non-singular big bang and big-crunch surfaces. The 
i'-dimensional causal structure is shown in the top right panel of Fig. [5) 

We have confirmed the existence of such solutions numerically. Solutions for p = 5 are shown in Fig. [71 where we 
have set Q = 1 and = M^j^^^dI"^- expected, there are solutions with 4 (top panel) and 2 (bottom panel) 
nodes for initial conditions at increasingly large (j). 
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FIG. 8: The extremal solutions connecting (j)- and 4> ^ oo (top panel) and cj>- and — oo (bottom panel). In the bottom 
panel, there is a singularity at finite r when a reaches a second zero, denoted by the asterisk. Near r —oo, the field 
approaches the potential maximum, and a approaches an exponential. Event horizons of the D-dimensional geometry are 
located at r = — cxd, where the two types of solutions in the top and bottom panels (going to </!> = ±cxd, i.e. R — or R — oo) are 
sewn together. For p odd, two regions of asymptotically flat D-dimensional space are matched across the horizon (two copies 
of the solution in the top panel), while for p even, a region of asymptotically flat D-dimensional space is matched across the 
event horizon to a region containing a timelike singularity (corresponding to the top and bottom panel). The causal structure 
is shown in Fig. [S] These are the non-dilatonic extremal p-branes of Ref. p^ . and the infinite range in r corresponds to the 
infinite proper distance separating any point from the event horizon. 

B. Flat FRW metric ansatz 

We now discuss the extremal solution, which is obtained for a flat FRW metric ansatz. Because the potential is 
negative in most regions of interest, the field cannot have a stationary point (/> = in a region with timelike r (from 
Eq. [3D] this would cause the Hubble constant to become imaginary). This, together with the condition a < from 
Eq. [3T1 implies that such spacetimes always begin with a singular big-bang and end in a singular big-crunch. If instead 
T is a spacelike variable, the evolution is in the upside-down potential —V, and it is possible to have a non-singular 
stationary point (j) = 0. As described in Sec. IIII B[ nonsingular solutions must be potential dominated as a — > and 
the metric coefficient a must behave asymptotically like a — exp(iJr) as r — > — oo, with 0(t — )■ — oo) = 0. The only 
two solutions that can meet these requirements interpolate between the inverted-potential maximum and ±oo, 
as shown in Fig. [8] These "extremal" solutions are the only trajectories where can loiter for an infinite range in r 
near a stationary point. All other solutions have a naked timelike singularity (the causal structure resembles the top 
left panel of Fig. [51). 

The non-singular surface at t = — oo corresponds to an event horizon in the _D-dimensional geometry. The infinite 
range in r between any point 4> and the potential maximum corresponds to an inflnite proper distance between any 
point in the D-dimensional geometry and the event horizon. The individual solutions interpolating between the 
potential maximum and (j) = ±c>o can be matched at the event horizon, corresponding to sewing together regions 
of asymptotically flat D-dimensional space and regions containing a timelike singularity. These full solutions are 
none other than the non-dilatonic extremal p-branes described in (15i |. These authors found that when p is odd, two 
identical regions of asymptotically flat space are connected across the event horizon, while for p even, a region of 
asymptotically flat space is connected to a region containing a timelike singularity across the event horizon. 

In order to determine which of our solutions can be pasted together, we can take the limit of the open FRW 
solutions as the endpoint of the motion approaches the inverted-potential maximum. At the stationary point in the 
region where r is spacelike, it is possible to analytically continue to a region where r becomes timelike, and the 
motion is about the potential minimum. For motion very near the minimum, we found above that the solution for 
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(j) was approximately a Gegenbauer polynomial with p + 1 nodes. For any finite displacement from the potential 
minimum, non-linearities cause the solution to become singular near the second stationary point in a, preventing a 
second analytic continuation back to a region where r is spacelike. For p even, the singularity occurs to the left of the 
minimum (0 has an odd number of nodes), and for p odd, to the right of the minimum (since <j) has an even number 
of nodes). Taking the limit where the stationary point of the field reaches the potential maximum, from the spacelike 
side, the fact that we started with an open universe is irrelevant since an infinite amount of proper time elapses before 
leaving the potential maximum. However, the properties of the timelike solution survive, and determine if the metric 
is invariant under reflection of about the maximum. Thus, for p odd, the spacetimes are completely non-singular, 
connecting two identical asymptotically flat regions of the spacetime across an event horizon. For p even, there is a 
timelike singularity behind the event horizon. This agrees with the analysis of Ref. [isj . 

C. Closed FRW metric ansatz 

Finally, we can briefly visit the closed FRW ansatz. Just as for the flat case, solutions with t timelike always 
begin with a singular big-crunch and end in a singular big-bang. For spacelike r the discussion in Sec. IIII 51 indicated 
that there is always a singularity at a = (which in this case will be a naked timelike singularity). This singularity 
need not be present if there are bouncing solutions where a is positive definite. From Eq. I31[ this requires that the 
field remain potential dominated over its entire evolution, which simply does not occur in the presence of the radion 
potentials (see Ref. Q for further discussion). We therefore conclude that all solutions with a closed FRW ansatz are 
singular. 

In conclusion, each of the solutions that we have found correspond to extended objects embedded in asymptotically 
flat D-dimensional space. The singular and oscillatory solutions are new classes of solutions to higher dimensional 
Einstein gravity, whereas the extremal solutions for p even or odd are identifled as the non-dilatonic p-brane solutions 
of Ref. [l^. We now turn to the case with a positive IJ-dimensional cosmological constant, A > 0, where again we 
will find new classes of solutions. 



V. SOLUTIONS FOR A > 

We now construct solutions where there is a positive cosmological constant A in the D-dimensional theory. The 
methods used are identical to those in the previous section, but the character of the solutions will be very different 
due to a number of qualitative changes in the radion potential. First, for any non-zero A, over some range of charge, 
the radion potential possesses two extrema (as outlined in Sec. Ill B[) . One is a maximum located at = and the 
other a minimum at = . Unlike the case where A = 0, the minimum can be either positive, zero, or negative 
energy depending on the charge Q. In addition, the potential in the case where A > approaches zero from above as 

(j) oo. 

The extrema of the potential give rise to AdSp+2 x S"^, A4p+2 x S"^, and dSp+2 x S"^ compactification solutions. 
There are also three qualitatively different classes of solutions that contain an asymptotically D-dimensional de Sitter 
region: the singular, interpolating, and Nariai solutions. The singular and interpolating solutions are generated from 
an open FRW metric ansatz, and the Nariai solution from a fiat FRW ansatz. There are important qualitative 
differences between the solutions in the case where the minimum of the radion potential is negative as opposed to 
zero or positive, and so we treat each of these cases separately. 

The possible trajectories when the radion potential has a negative minimum are shown in Fig. [51 the possible 
solutions for a positive or zero minimum are shown in Fig. |5l and the causal structure of the corresponding D- 
dimensional spacetimes is shown in Figs. [TI] and 1121 Notable among these geometries are the interpolating and 
Nariai solutions, which have the remarkable property of interpolating between Z?-dimensional and effectively p + 2- 
dimensional regions across event horizons. We summarize the cosmological significance of these solutions at the end 
of this section. 

For fixed charge, there is a compactification solution for each of the critical points of the potential (0, 1, or 2 
critical points are possible). Positive, zero, and negative energy minima of the radion potential yield a dSp+2 x S**, 
Mp^2 X S"^, and AdSp+2 x S** geometry respectively. Since these solutions sit at a potential minimum, they are stable 
to homogenous perturbations in the radion field. The compactification solution produced by sitting at the radion 
potential maximum will always be dSp+2 x S'^ . This solution is unstable to homogenous perturbations in the radion 
field, which would cause spacetime regions to evolve away from the maximum at 0+ either towards smaller or larger 

The construction of the Anti-de Sitter compactification solutions will proceed as described above. For positive or 
zero energy minima, the behavior of the scale factor is different. Minkowski space can be foliated into flat or open 
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FIG. 9; A sketch of the solutions for A > when the radion potential has a negative energy minimum. Motion above the 
potential occurs in a region where r is a timelike variable, and motion below in regions where r is a spacelike variable. An 
asterisk denotes the occurrence of a curvature singularity along the trajectory. Non-singular stationary points are denoted by 
the filled circles. All of the solutions contain a singularity except for the compactification solutions that sit at the critical points 
of the potential for all r. At the minimum, an AdSp+2 x geometry results, and at the maximum a dSp+2 x S"' geometry 
results. For an open FRW ansatz, the interpolating solution has two non-singular stationary points and connects <^ ^ cx3 to 
a region containing a spacelike singularity. All other choices for the location of the stationary points will lead to a timelike 
singularity. In either of these cases, the D- dimensional solution approaches de Sitter as r (x>. For a flat FRW ansatz, the 
Nariai solution connects regions interpolating between 0+ and (j) — > ±oo over an infinite range in r. The Nariai solution has a 
stationary point at 0+ , but the geometry is distinct from the compactification solution (although the near-horizon geometry of 
the Nariai solution is locally dSp+2 x S''). 



FIG. 10: A sketch of the solutions for A > when the radion potential has a zero or positive energy minimum. There is 
one de Sitter and one de Sitter or Minkowski compactification solution at the critical points of the potential. As for the case 
where the radion potential has a negative minimum, for an open FRW ansatz there is an interpolating solution that has two 
non-singular stationary points connecting the asymptotic — > oo region to the basin of attraction of the potential minimum. 
Motion away from the second stationary point will lead to a completely non-singular evolution towards </>_, asymptotically 
producing a region of p + 2 dimensional Minkowski or de Sitter space. Solutions with any other stationary points contain a 
timelike singularity. For a fiat FRW ansatz, the Nariai solution connects regions interpolating between and <;/>_ or <j!> ~* oo 
over an infinite range in r. The Nariai solution is in this case completely non-singular. 

spatial sections, and de Sitter space flat, open, or closed spatial sections. In any case, the equation of motion for the 
scale factor in regions where r is timelike is again given by 
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FIG. 11: Causal diagrams for A > solutions with an asymptotically D- dimensional region and a radion potential with a 
negative minimum. Surfaces of constant r are shown as the light solid lines. Starting from the left, we show the singular, 
interpolating, and Nariai solutions. All solutions contain a timelike or spacelike singularity separated from asymptotically 
D- dimensional de Sitter space (denoted by the D) by event horizons. The horizons are in each case color coded with the 
stationary points of the radion field trajectories in Fig. [O] We highlight the solutions relevant for dynamical compactification 
in the shaded box. 



where 



^a^YJp f , -IV 61) 

The initial conditions are set by the Friedmann equation Eq. [3D1 and depend on the spatial curvature. 



A. Open FRW metric ansatz 



Beginning our study of dynamical solutions with an open FRW metric ansatz, we need to consider the evolution 
of the field in three separate regions of the potential: the asymptotic (/> ^ oo region, the region under the potential 
barrier, and the region in the basin of attraction of the minimum. For A > the potential approaches zero from above 
as — *■ oo, implying that r is a timelike variable in the asymptotic region. If there is evolution under the barrier (this 
requires that the energy density in the field is small enough), r must be a spacelike variable. Finally, for evolution in 
the basin of attraction of the minimum, r is again a timelike variable. The evolution in each of these three regions 
can be sewn together across non-singular big-bang and big-crunch surfaces where = and a — 0, just as in the 
previous section. 

In the asymptotic region (p — > oo, where r is timelike, we can again look for attractor solutions that are approached 
starting from generic initial conditions. The dominant term in the potential at large is 



m - 4r exp ( J ] . (62) 



Ml \ y Pip + q) Mp+2 

Again, wc will find it convenient to work in terms of the dimensionless variables defined in Eq. 1351 From the asymptotic 
form of the potential, we can identify the constants 



A j q 4 ^ p + q q 2A 



^ M2' ^ ^\lp{p + q)' [pi^) ^ q p\l Ml{p + q){p + q+iy ^^^^ 

In this case, wc have (5 > 1, and the field dynamics dominate the contribution from curvature. The metric in the 
i'-dimensional Einstein frame is given by 

Mlds' = (cf)-2 r_^-2 ^ -2(^-)2(p+g)/g ^^^2 ^ ^^^^^2 ^^^2^1 ^ (^f )2p/«df]2_ (g^) 
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FIG. 12: Causal diagrams for A > solutions with an asymptotically D-dimensional region and a radion potential that has 
a positive or zero minimum. Surfaces of constant r are shown as the light solid lines. Starting from the left, we show the 
singular, interpolating, and Nariai solutions. In the interpolating solution, future/past infinity is split between non-singular 
regions of different effective dimensionality (the dimensionality is labeled D or p + 2). When the minimum has zero energy, the 
p+ 2-dimensional parts of future infinity is null (denoted by the dashed "hat"). The Nariai solution connects the contracting 
portion of a 7?-dimensional de Sitter space to a non-singular p -|- 2-dimensional flat FRW universe across an event horizon. The 
horizons are in each case color coded with the stationary points of the radion field in Fig. (5] We highlight the solutions relevant 
for dynamical compactiflcation in the shaded box. 



We can define a new time coordinate 

cMdT ^ log(cf), 

in terms of which (after substituting for the appropriate values of c from Eq. I47p the metric becomes 



flg CXP 




sinh^ Xdnl + dfll] 



(65) 



(66) 



This is the metric ofD=p + q + 2 dimensional de Sitter space with cosmological constant A (this analysis was first 
performed in Ref. Q). 

As in the previous section, if a(T = 0) = at a stationary point 0(t = 0) = 0, then there is a corresponding event 
horizon in the asymptotically i?-dimensional region. Otherwise, there is a naked spacelike singularity corresponding 
to a singular big-bang or big-crunch surface. In the non-singular case, we can continue across the stationary point to a 
region where the field evolves in the inverted potential and r is spacelike. Considering motion about the near-vicinity 
of the inverted maximum, the scale factor will have a sinusoidal behavior, and the field equation has a solution in 
terms of Gegenbauer polynomials (as was the case for motion about the potential minimum for A = 0). The number 
of nodes is given by the positive root of 



cr^ + (p + l)cr = 



p{p+l)\V"{^+)\ 



2 |V^(0+)l 

In the limit of small Q, we can use Eqs. [55] and [25] to find the ratio 

l^"(0+)l _ 4 
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Solving for a, 
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FIG. 13: The interpolating solution for p = 2, q = 2, Q — 3.2, and A/M^ — .1 generated from the initial condition (j)o — 2.01. 
The field interpolates between either side of the potential barrier while the scale factor evolves between its two zeros. 



This ranges between 1 < a < 2 for all p and is never integer. 

Starting from an arbitrary stationary point to the right of the potential maximum yields a singular solution, since 
a is not integer. The Z?-dimensional causal structure of such solutions is shown in the left panel of Figs. [TT] and [T2l 
However, this singularity is potentially avoidable since the field has a sufficiently long time to cross the potential 
maximum at least once before a reaches its second zero. Increasing the initial displacement from the maximum, 
the frequency of oscillations in the field decreases, and there will always be a critical point at which the traversal 
time of matches the period of the scale factor. This is essentially Coleman's undershoot/overshoot argument for 
finding instantons in a scalar field potential [25l |. For different values of the charge, there is a unique stationary point 
that yields such a non-singular "interpolating" solution. Although the requisite set of initial conditions seems rather 
special, we will find in Sec. IVII that there are instanton solutions that pick them out uniquely. 

We have confirmed the existence of such solutions numerically. An example for the parameters p — 2, q — 2, 
Q = 3.2, and K/M'jj = .1 is shown in Fig. [131 Here, it can be seen that the radion field interpolates between either 
side of the barrier as the scale factor evolves between its two zeroes. In the Z?-dimensional geometry, this describes 
the behavior of the g-sphere between two event horizons. As the value of the charge increases for fixed p, q, and A, 
the potential maximum broadens (recall that the maximum and minimum of the potential merge at a finite charge 
Eq. [M]). At a large enough value for the charge, we have cr < 1, implying that the field does not have sufficient time 
to cross the barrier before the scale factor reaches its second zero. In these cases, the interpolating solution docs not 
exist. 

At the second non-singular stationary point {(p < 4>+) of the interpolating solution, we can continue to a region 
where r is again a timelike variable. The properties of the solution in this region depend on the sign of the potential 
minimum. For small charge (much less than the the limit Eq. I24p , the potential near the minimum is negative and 
approximately the same as for A = 0. Evolution away from the stationary point of the interpolating solution always 
results in a singularity as the non-zero kinetic energy of the field blue-shifts in the negative vacuum energy background 
(this is in exact analogy with the big-crunch that occurs inside of negative energy vacuum bubbles [26|). If instead 
V{(j)-) > 0, the field will simply experience damped oscillation around the positive minimum while the scale factor 
grows monotonically. The solution in this case is completely non-singular. 

Solutions of this type interpolate between the asymptotic Z?-dimensional region (f> oo and the basin of attraction 
of the radion potential minimum where the spacetime is effectively p + 2-dimensional. The causal structure of these 
solutions is shown in the center panel of Figs. [11] (for a negative minimum) and 1 121 (for a zero or positive minimum). In 
each diagram, there are three qualitatively different regions corresponding to the three segments of the field evolution 
described above. In each case, there is a patch in which the field evolves io (p ~* oo, asymptotically reaching D- 
dimensional de Sitter space (the right portions of past and future infinity marked "D" on the causal diagrams) and a 
region where the field interpolates between non-singular stationary points under the potential maximum (the central 
diamond on the causal diagrams) . There is then a region that evolves in the basin of attraction of the radion potential, 
ending in a singularity if the minimum is of negative energy (as in Fig. Ilip or in a p -f- 2-dimensional asymptotically 
fiat or de Sitter region (the left portions of past and future infinity marked "p-f2") if the minimum is zero or positive 
energy (as in Fig. I12p. At the boundary of each segment, there is an event horizon. 

The regions evolving in the basin of attraction of the radion potential are p + 2-dimensional open FRW universes 
with a non-singular big-bang surface. This is extremely interesting from the standpoint of cosmology, since there 
exists a pre-big bang epoch in which the effective dimensionality of spacetime changes. 

What happens to the interpolating solutions when the charge goes to zero? In this limit, the minimum of the 
potential goes to — s- — oo. At Q = 0, the minimum disappears, since the positive contribution to the potential 
Eq. [201 is no longer present. Studying a few numerically generated examples, we find that the stationary point to 



21 



the left of the potential maximum follows the location of the minimum towards (f> — > — oo. At Q = 0, the location 
of the stationary point will be at i? = 0, which is just the non-singular origin of coordinates. The spacetime has a 
single horizon, corresponding to the stationary point to the right of the maximum. This is nothing other than empty 
ZJ-dimensional dc Sitter space. 

Finally, we note that in regions with timelike t it is possible to find a portion of the oscillatory solutions described 
in the previous section. For small charge, the potential in the vicinity of the minimum is largely unaffected by the 
positive cosmological constant, and so the analysis will proceed exactly as before. As the charge increases, the depth 
of the minimum is decreased, and the term in the potential proportional to A begins to contribute to the second 
derivative. This causes the ratio Eq. [58] to diverge as V{(f>^) — > 0. As this occurs, it becomes possible to find an 
increasing number of oscillatory solutions since the frequency of the scale factor oscillations goes to zero while the 
frequency of field oscillations stays roughly constant. The number of solutions diverges as ^ — *■ 0. However, the non- 
singular stationary point of the interpolating solution never match the non-singular stationary point of the oscillatory 
solution. This means, in contrast to the solutions with A = 0, that it is impossible to match these oscillatory solutions 
to an asymptotically D-dimensional region. 

B. Flat FRW metric ansatz 

For a flat FRW ansatz, the extremal solution evolves away from the potential minimum in a region where r is 
spacelike, and can only exist when V{(f>-) < 0. However, unlike the case where A = 0, these solutions cannot evolve 
to — > CX3 (recall that the interpolating solution is the only non-singular solution in the region under the potential 
barrier). Instead, these solutions have two timelike singularities separated by an event horizon. 

In the case of the Nariai solution, the field evolves away from the potential maximum towards either 0_ or — > oo 
over an infinite range in r. It is similar to the extremal solutions found for A = 0, except that here r is everywhere 
timelike. A numerically generated example for a radion potential with positive minimum is shown in Fig. 1141 This 
solution is completely non-singular, and corresponds to a contracting D-dimensional de Sitter space separated by an 
event horizon from an asymptotically p + 2-dimensional de Sitter space with a g-sphere of constant size. Unlike the 
extremal solutions of the previous section, an asymptotically i?-dimensional region where </> ^ oo is always matched 
across the horizon to an asymptotically p + 2-dimensional region where — *■ 0_ (if the radion potential is zero or 
positive) or a singularity at — > — oo (if the radion potential minimum is negative). This is due to the fact that the 
Gegenbauer index a for small- amplitude solutions is always less than 2, and therefore motion across the potential 
maximum never has two nodes. 

Like the interpolating solutions, the Nariai solutions interpolate between an asymptotically D-dimensional de Sitter 
space and a p + 2-dimensional region that is either singular (when the minimum of the radion potential is negative) 
or asymptotically de Sitter/Minkowksi (when the minimum is positive or zero energy). The causal structure is shown 
in the right panel of Figs. [TT] and fT2l Unlike the interpolating solutions, future and past infinity are not bifurcated 
into separate regions. 

C. Summary of the cosmologically interesting solutions 

The compactification, interpolating, and Nariai solutions all have spacetime regions that contain ap-|-2-dimcnsional 
FRW universe with a stabilized g-sphere. This allows for the intriguing possibility that our observable universe could 
actually be embedded in one of these geometries. 

Compactification solutions generated at the minimum of the radion potential are equivalent to one of the vacuum 
geometries AdSp+2 x -S', Aip+2 x S'^, or dSp+2 x -S"'. These solutions are stable vacua, and the existence of de 
Sitter solutions is compatible with the observed cosmological acceleration in our universe. The compactification 
solution generated at the radion potential maximum is also a dSp+2 x 5"' geometry, but it is in unstable equilibrium. 
Perturbations will cause the radion field to evolve away from the potential maximum in different spacetime regions. 
Some such regions can evolve towards the minimum of the radion potential, resulting (at least locally) at late times 
in a p -|- 2-dimcnsional de Sitter vacuum. 

The Nariai solutions (which are generated from a flat FRW metric ansatz) interpolate in a timelike manner from 
a contracting D-dimensional de Sitter space to an expanding infinite p + 2-dimensional flat FRW universe that can 
evolve towards zero or positive vacuum energy. These regions are again separated by an event horizon, although in 
this case it extends across the entire universe. The radion field is located at the potential maximum on the event 
horizon, and so the near-horizon geometry will be locally equivalent to the compactification solution dSp-1-2 x S*"*. 
Since the radion field is dynamical, there will also be portions of the Nariai solution that arc locally equivalent to the 
spacetime that results from perturbing the compactification solution away from the potential maximum. 
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FIG. 14: An example ol the Nariai solution for p = 2, q — 2, Q — 3.2, and A/M|, = .1. The field evolves from the potential 
maximum to either (f)- (top) or ^ oo (bottom) while the scale factor a monotonically increases. The two solutions are 
matched across the event horizon of the D-dimensional geometry located ai t — —oo. This connects the contracting portion of 
a D- dimensional de Sitter space to a flat p + 2-dimensional FRW universe that asymptotically has positive vacuum energy. 



The interpolating solutions (which are generated from an open FRW metric ansatz) contain an infinite p -\- 2- 
dimensional open FRW universe that can evolve towards positive vacuum energy at late times. The big bang of the 
FRW universe is non-singular, and corresponds to an event horizon in the D-dimensional geometry. On the other 
side of the event horizon, there is a region in which the radius of the g-sphere grows in a spacelike manner. As the 
radius of the g-sphere grows, another horizon is eventually reached, across which is an asymptotically _D-dimensional 
de Sitter region. 



VI. DYNAMICAL COMPACTIFICATION 



Thus far, we have discussed the classical solutions of Z?-dimensional Einstein gravity with a cosmological constant 
and flux. We have identified a number of solutions that contain a lower-dimensional cosmological spacetime, but these 
are full time-symmetric solutions, and do not represent a transition. We now examine the semi-classical properties of 
the theory, where transitions to these classical solutions can occur from empty de Sitter space. 

This process has an analog in black hole physics. Because of its thermal character, it has long been appreciated 
that four-dimensional de Sitter space is unstable to the nucleation of both charged and uncharged black holes. Such 
nucleation events can be described in the semi-classical approximation by Euclidean instantons, with the probability 
per unit four- volume given by 

F = Ae-('^'"=*-^<^«), (70) 

where A is a prefactor, Smst is the Euclidean action of the instanton, and Sds is the action of Euclideanized de Sitter 
space. In higher dimensional de Sitter space with g-form field strength, given the proliferation of classical solutions 
discussed in previous sections, one might expect that there is a finite probability for their nucleation. 

Indeed, we will find that this is the case. Examining the solutions for A > presented in Sec. |Vl we must identify 
those solutions that when Euclideanized yield a finite action. A sufficient condition for finite action is that the 
Euclidean manifold is compact, and so the most obvious candidates are the dSp+2 x S"^ compactification solutions. 
Euclidean de Sitter space of dimension p -I- 2 is ap-|-2-sphere, which yields the compact Euclidean manifold 5^+^ x S*^. 
Moving to the solutions with an open metric ansatz, we note that they possess an 5*0(^-1-1, 1) symmetry. The obvious 
analytic continuation is to a metric with 0{p + 2) symmetry, which can be related to the metric for spacelike r Eg. 1331 
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FIG. 15: The spacetime picture of dynamical compactification when it is mediated by the Euclidean interpolating (left) and 
compactification (right) solutions. On the left, the initial D- dimensional de Sitter space is matched across the nucleation surface 
to the interpolating solution of Sec.|Vl The surfaces of constant <j!> are shown as solid lines, and interpolate across the barrier of 
the radion potential along the nucleation surface. On the right, roughly a horizon-volume of the initial D-dimensional de Sitter 
space fluctuates to the compactification solution at the maximum of the radion potential. After the fluctuation, the radion falls 
from its unstable equilibrium in different spacetime regions, producing pockets where the field either evolves to the minimum 
of the radion potential or off to infinity where the g-sphere decompactifles. 



by taking x ^ «(x + 37r/2), yielding 



dr^ + a(r 



(71) 



This manifold is also generally compact when the scale factor a(T) is bounded. Because the radial coordinate x does 
not appear in the equations of motion for a and cj) (Eqns. [501 and 1^5)1 . they remain unaltered. The Euclidean profiles 
for a and (p are then identical to those presented in Sec. |V] for regions with spacelike r (the region under the potential 
barrier). Therefore, another class of finite action instantons can be obtained from the Euclideanized interpolating 
solutions in the region underneath the potential barrier, where air) possesses two zeroes (an example is the profiles 
shown in Fig. [T3|) . Finally, we can examine the Nariai solution for a flat FRW metric ansatz. In this case, there is no 
analytic continuation to a compact manifold, and so no finite action action instantons will result. 

The finite action instantons will mediate a transition from empty D-dimensional de Sitter space to either the 
compactification or interpolating spacetimes. The spacetime picture of these two processes are shown in Fig. 1151 
Because in each case there is a region after the transition that relaxes to a p + 2-dimensional vacuum, this can rightly 
be viewed as a dynamical mechanism for the compactification of q dimensions. 

When dynamical compactification is mediated by the Euclideanized interpolating solution, the spacetime resembles 
the left panel of Fig. 1151 A /^-dimensional de Sitter space is matched onto the Lorcntzian interpolating solution 
across a nucleation surface, which by time translation, can be placed at the throat of the background de Sitter space. 
Evolving, there are two infinite open FRW universes separated by event horizons, one of which is p -t- 2-dimensional. 
This process splits future infinity between a p -I- 2 and D dimensional region. 

If the Euclidean compactification solution is the relevant instanton, then the spacetime picture is sketched in the 
right panel of Fig. [151 The interpretation is very similar to that of the Hawking-Moss instanton for a scalar field 
coupled to gravity [27], or the fiuctuation of a Nariai black hole Q. In particular, a roughly horizon-sized patch of 
the D-dimensional de Sitter space fiuctuates into a solution very close to the Lorentzian compactification solution at 
the maximum of the radion potential. However, this solution is unstable, and will fall either into the asymptotically 



^ It is also possible to obtain a finite Euclidean action from the compactification solution at a positive minimum of the radion potential. 
However, since there are no negative modes, this cannot be interpreted as mediating a transition. For a discussion of this point, see e.g. 
Ref. m. 
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D-dimensional region at (/) ^ oo or towards the p + 2-dimensional vacuum at (/<_. This evolution will locally be similar 
to the Nariai solution described in the previous section, although the global causal structure will be different. 

In addition to dynamical compactification, the inverse process can also occur, whereby ap + 2 dimensional vacuum 
can undergo a "decompactification transition" to a region that empties out into _D-dimensions [2I, [2^ . This is medi- 
ated by the same instantons responsible for dynamical compactification [2:]. The spacetime picture for spontaneous 
decompactification is similar to Fig. [151 with the region of Z?-dimensional de Sitter space to the past of the nucleation 
surfaces replaced by a p + 2-dimensional de Sitter. 

In the next subsection, we study the nucleation rates in a general radion potential. We then briefly describe 
the global spacetime picture that emerges from considering multiple transitions between a U-dimensional de Sitter 
space and the great variety of possible p + 2-dimensional vacua. We largely avoid interpreting our results in terms of 
statistical predictions for physical quantities, but simply present the possibilities and in some cases relative likelihoods 
of populating different vacua. 



A. Populating a landscape of vacua 

The many radion potentials for fixed p, q, and A that exist for different values of the charge define a landscape of 
lower dimensional vacua. Each of the potential minima have a different p + 2-dimensional vacuum energy, ranging 
from arbitrarily negative values at Q — > either to zero as Q — > 00 in the case where A = or to a bounded positive 
value as Q — > Qmax (by Eq. I24p in the case where A > 0. If the charge is quantized, then there are an infinite number 
of vacua for A = where Q can become arbitrarily large, but only a finite number for A > where the charge must 
be less than Qmax- Allowing for the existence of multiple types of g-form flux by replacing 

we can also flnd collections of vacua with different effective dimensionality. The dynamical compactiflcation mechanism 
described above provides a means to evolve from D-dimensional de Sitter space into each of these vacua. In this section, 
we will quantitatively compare the nucleation rates for each transition. 

In analogy with the nucleation rate for black holes in four dimensions Eq. 1701 we take the probability per unit 
£)- volume of the background de Sitter space to nucleate a solution to be 



r = A cxp 



-{Sinst - 5*^5^) • (73) 



The prefactor A is of mass dimension D and represents the first quantum corrections to the probability. This quantity 
is in general difficult to calculate (requiring us to step away from the minisuperspace of metrics we have considered), 
but we expect it to be of order A^/^. The action for the higher dimensional background dS space is given by 

^(D) ^ _ / d^x.raiR - 2A) = -(v + a+ l)(P+'^+^y^h + a)^P+'^y'Vom,.^„^.) ( — 



j d''x^{R-2A) = -ip+q + iYp+''+'y\p+qYp+'^y'Yo\{n,+,+2)i^^j ■ (74) 



Substituting the metric ansatz Eq. [7T] for the instanton into the action Eq. [11] and integrating over the angles of the 
p + 1-sphere, the Euclidean action is given by 



Sinst = -2Yo\{np+i) J dra"^ 



y I 1 P(P+1) 
2 a2 



(75) 



There are two types of instantons that can be responsible for dynamical compactification: the Euclideanized 
compactification and interpolating solutions. Both types will have an action given by Eq. 1751 It is possible to find the 
compactification solutions analytically, but it will be necessary to construct the interpolating solutions numerically, 
and we will therefore find it convenient to work with the dimensionless variables defined in Eq. 1351 A characteristic 
scale for the potential is LUa, which is essentially the barrier height, and is defined in Eq. l61l Substituting for the value 
of V{4>+), LOa is given by 

Vol {nq f^" [p +l){p + qY'P+^[q - 1)9/P V ' 
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Setting ti = M, 



P/(P+2), 2/(p+2) 
p+2 



, the dimensionless FRW variables are given by 

V 



V = 



p+2 



T = WaT, a = (jJatt, 



(77) 



Substituting into the action, we obtain 

p 



Si', 



-2Vol(r!p+i) 



Mp+2 



V - 



1 p{p + 1) 



2 a2 



-2Vol (f^p+i ) Vol (17,) + 1)^/2 ( ^ 



2 \ (P+'?)/2 



V - 



1 p{p + 1) 



2 a2 



.(78) 



We use this form of the rescaled action throughout the rest of this section in our computation of the instanton action. 
A convenient parametrization of the nucleation rates is in terms of the ratio 



^ - C IqiD) 

" — i-'mst/ "J • 
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The nucleation probabilities are then given by 

r = ^ exp 

which grows with increasing a since 

'S'ds^ < 0- For the semi-classical formalism to accurately describe the transition 
process, we must always have suppressed transitions, requiring a < 1. We will find that this is indeed the case for 
all of the solutions we construct. Unless a ~ 1, the rates are set by 5*^5^ which is negative and large in magnitude 
unless A ^ M^. For such a large A, we expect the semi-classical description of the background D-dimensional de 
Sitter space to break down. 

Substituting the instanton action Eq. [78] and the background de Sitter action Eq. [Til we find that 



Vol (r!p+i) Vol (Qg) {p + l)f/2(q - l)'?/2 



voi(r2j 



'P+q+2 



) (p + g+l)(P+9+2)/2 



- l p(P+l) 

^2 a2 



(81) 



The first (and simplest) finite-action instanton we consider is the Euclideanized compactification solution, where 
the field sits at the maximum over the entire evolution in r. The solution for the scale factor is air) = ui~^ sin(a;oT). 
Substituting the dimensionless time variable f = cJaT, dimensionless scale factor a(f) ~ sinf, and dimensionless 
potential V = V/{Mp^2P^a) — ~P{P + into the instanton action, we obtain 



S^nst = -p{p + l)Vol [VLp+i) \-^] J dfsmP-' f [1 - sin2 f] 



The integral can be evaluated in terms of F functions, and after some manipulation we find 



S^„st = -(p+ l)V0l(f]p+2) 



(82) 



(83) 



This is equivalent to the volume of a p -I- 2-sphere of radius u!~^, as expected from the fact that we are simply 
Euclideanizing a, p + 2-dimensional de Sitter space. For small Q, we can substitute for Ua using Eq. 1761 The value of 
a in this case is 



Vol (»p+2) Vol (»,) {p + l)P/2+l(g - 1)^/2 

voi(r!p+,+2) (p + q + i)(p+'?+2)/2 ■ 



(84) 



In Fig. [TBI we plot the value of a in Eq. [84] (valid for small Q) for fixed D over a range of p. Interestingly, this 
function peaks at a value of p specified by the total dimensionality. This is due to the relative volumes of the unit 
spheres. Relevant for our universe is that for D = 8, the nucleation rate for Euclidean compactification solutions with 
p — 2 (and therefore four non-compact dimensions) is largest. There are also marginal cases: for D = 7, the rate for 
p = 1 and p = 2 are identical while for D = 9, the rate for p = 2 and p — i are identical. Also of note is the fact that 
a is in general larger when the total dimensionality is increased. 

For appreciable values of Q, the relative importance of the various terms in the potential changes, and Ea. l26l becomes 
a poor estimate for V{(j)+). In Fig. [T7] we plot the exact value of a versus Q (circles) for {p = 2,g = 2, A = .1} 
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FIG. 16: The value ol ct (Ea. l84|l . the ratio of the instanton action to the D-dimensional de Sitter action, for fixed D and various 
p in the limit of small Q. Larger values of a correspond to a faster nucleation rate. Circles, squares, diamonds, triangles, and 
inverted triangles are for D = 10, 9, 8, 7, 6 respectively. In each case, p can range from to 1 to D — 2. The value of a peaks at 
a value of p dependent on the total dimensionality. 

determined by numerically finding V^(0+). Here, it can be seen that at small Q, a approaches the value determined by 
Eq. [M] (depicted by the solid horizontal line). Further, a decreases with increasing Q. For Q larger than the critical 
value defined in Eq. 1241 there are no potential extrema, and the Euclidean compactification solutions do not exist. 

Now, we consider the second instanton, found by Euclideanizing the interpolating solution. Evolution in the presence 
of the Euclidean potential was already discussed in Sec. |Vl This generated the portion of the interpolating solution in 
regions with spacelike r (underneath the barrier). The profiles for a(r) and 4>{t) are qualitatively identical to those 
depicted in Fig.[T31 where we set (/) = at the endpoints in a, and the field interpolates across the potential barrier. 
It is possible to construct such solutions numerically, and compute the value of a using Eq. 1811 

In Fig.[T71 we plot the value of a versus Q for the Euclidean interpolating solutions (squares) using {p = 2, g = 2, A = 
.1}. The values of a at any fixed Q are always larger for the Euclidean interpolating solutions than for the Euclidean 
compactification solutions, implying that the nucleation rate is also larger for the Euclidean interpolating solutions. 
As Q increases, a decreases, and the action of the interpolating solutions approaches that of the compactification 
solutions. At small Q, a is very close to unity. Because the potential minimum goes to V- — oo as Q ^ 
(see Eq. [231), our numerical solutions become unreliable in this regime, and it is not possible to determine what the 
asymptotic value at Q = is. However, as we discussed in SeclVj the zero charge limit of the interpolating solutions 
should be empty D-dimensional de Sitter space, which would be consistent with a — > 1. At Q = 0, the beginning and 
end states of the "tunneling" process are identical. 

The properties of the potential maximum are altered as Q increases. Eventually the index of the Gegenbauer 
polynomials defined in Eq. [57] becomes less than one, and it is impossible to find an interpolating solution. This 
generally occurs for Q less than the value for which there are no critical points of the potential (defined in Eq. IM)) , 
and thus there is a range in Q over which the Euclidean compactification solution is the only viable transition 
mechanism. This is depicted by the dark (green) shaded region in Fig. 1171 In general (for arbitrary p, q, and A), 
there will be a window of values of Q with a positive cosmological constant minimum and a Euclidean interpolating 
solution, depicted by the light (blue) shaded region of Fig. [TT] 

In Fig. [181 we zoom in on the shaded region of Fig.[T7l As Q is increased, the endpoints of the interpolating solution 
approach the potential maximum, eventually merging when it becomes impossible to find the interpolating solution. 
The action for the interpolating solution approaches that of the compactification solution, dominating the decay rate 
when it exists. 

To the extent that our dimensional reduction captures the relevant physics, the number of negative modes for 
each type of solution can be determined by the arguments of Ref. [1^. The interpolating solution (similar to the 
Coleman-de Luccia instanton (soj). when it exists, will always have a single negative mode, and therefore represents a 
semi-classical instability. The Euclidean compactification solution (similar to the Hawking-Moss instanton [23|) will 
have a single negative mode only when the interpolating solution does not exist (and it therefore dominates the decay 
rate as illustrated in Figures [T71 and [T5| . 

It is also possible to compare the actions for interpolating solutions with variable p in a fixed total number of 
dimensions. In Fig. [191 we plot a vs the normalized charge Q/Qmax in D = 8 for p = 1, 2, 3, 4 (numerical data points 
are shown as circles, squares, diamonds, and triangles respectively). The value of Q yielding an approximately zero 
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FIG. 17; A plot ol a versus Q for the Euclidean interpolating (squares) and compactification (circles) solutions with {p = 2,q = 
2, A = .1}. The value of a for the compactification solutions approaches the value defined in Eg. 1841 (denoted by the horizontal 
solid line at a = .72) at small Q. A few important thresholds in Q are denoted by the shaded regions. For 3.162 < Q < 3.35, 
there exists a minimum of the radion potential that has positive energy. Between 3.162 < Q < 3.35, in the light (blue) shaded 
region, both the Euclidean interpolating and compactification solutions are allowed. For 3.35 < Q < 3.65, in the dark (green) 
shaded region, only the Euclidean compactification solution exists. When it exists, the interpolating solution always makes the 
dominant contribution to the decay rate. 
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FIG. 18: The shaded region of the q vs Q plot of Fig. 1171 in which the radion potential has a positive minimum and the Euclidean 
compactification (dashed line) and interpolating (solid line) solutions exist. As Q is increased, the instanton endpoints approach 
the potential maximum, and the value of a for the interpolating solution approaches that for the compactification solutions 
from above. When it exists, the interpolating solution always makes the dominant contribution to the decay rate. 



energy p + 2 dimensional vacuum in each case is denoted by the large red circle. At fixed Q/Qmax, the solutions for 
p = 2 always have the largest value of a (and hence the highest nucleation rate). It is also of interest to compare the 
rates to zero energy p + 2 dimensional vacua, in which case we find that the p — 2 and p — 3 solutions have comparable 
values for a (because of the numerical error incurred in doing the computation it is unclear which dominates in this 
case). In general the comparison of nucleation rates between solutions is complicated, since there is a value of Q 
at each p for which the nucleation rate is the same, and therefore vacua of different dimensionality and different 
cosmological constant will be populated with identical frequency. It is necessary to determine some physical criterion 
for making the comparison, a topic we return to in the next section. 

The formation of the interpolating solution brings into being a region that relaxes to a p + 2-dimensional big crunch 
(when the vacuum energy is negative), or an open universe with zero or positive vacuum energy. However, when 
the radion potential has a positive minimum, the instantons we have found can also mediate a decompactification 
transition [3, [2^ . The late-time p + 2-dimensional de Sitter region produced in dynamical compactification will 
therefore in some patches undergo transitions back to ZJ-dimensional de Sitter space. The tunneling rate out of a 
metastable p + 2 dimensional de Sitter vacuum is given by 



r = Aexp -{S.. 



q(p+2)n 



) = A exp 



o(p+2) 



(1-/3) 



(85) 
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or 
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FIG. 19: The value of a in 7? = 8 for p — 1, 2, 3, 4 (blue circles, purple squares, yellow diamonds, and green triangles 
respectively) versus the normalized charge Q/Qmax- The red circle denotes the value of Q yielding an approximately zero 
energy p + 2 dimensional vacuum for each p. 
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FIG. 20: A plot of /3, the ratio of the instanton action to the p+ 2-dimensional de Sitter action, versus Q for the interpolating 
(solid line) and compactification (dashed line) solutions with {p = 2, g = 2, A = .1}. The values of P are nearly identical for 
the two solutions. At Q ~ 3.162, the minimum of the radion potential goes to zero, sending /3 0. The interpolating solution 
no longer exists for Q > 3.35, and at the upper range of the plot Q = 3.65, there are no potential extrema. 



where the action of the lower dimensional de Sitter vacuum is 
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and the quantity /3 (a ratio of actions analogous to a) is given by 
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A plot of P versus Q for 2, g = 2, A = .1} is shown in Fig. [20l As ¥{(/)-) 0, 5, 



instanton action approaches a constant in this limit, P 

c(P + 2) 

space is extremely suppressed, with F ^ e . In the limit of a vanishing p 
the vacuum is completely stable. 
In summary. 



(p+2) 
dS 
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— oo. Because the 
2-dimensional dS 



0, and the tunneling rate out of the p 

2-dimensional cosmological constant, 



• The value of a at a fixed normalized charge Q /Qmax is maximized for a particular p or range of p depending 
on the total dimensionality. If D — 7,8,9 then solutions with p = 2 either maximize a (for D = 8) or are 
comparable to solutions with p — 1 (for D = 7) or p = 3 (for D = 9). A detailed comparison of nucleation rates 
for vacua with different dimensionality and cosmological constant should be possible. 
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• For fixed p and q, the rate is highest for small Q, and therefore the lowest (negative) energy vacua. 

• When it exists, the rate is higher at fixed Q, p, and q for the interpolating solution than the compactification 
solution. There is in general a range of Q for which there is a positive energy minimum, but no interpolating 
solution. 

• Because the instanton action is finite as the vacuum energy of the minimum goes to zero, the rate into a p + 2- 
dimensional vacuum with a small (positive) vacuum energy is much higher than the rate out. p + 2-dimensional 
Minkowski vacua are completely stable. 



B. Global structure of the multiverse 



Given the existence of the instantons described in the previous section, empty D-dimensional de Sitter space in a 
theory with g-form field strengths will inevitably populate all of the lower-dimensional vacua of the theory. Because the 
probability per unit 13- volume to undergo dynamical compactificaton is much smaller than a D-dimensional Hubble 
volume H^^T <^ 1, it is very rare to have more than one transition per Hubble volume per Hubble time. This implies 
that the D-dimensional de Sitter space can never be totally eaten up by transitions and can populate many different 
lower-dimensional vacua, in exact analogy with four-dimensional eternal inflation. If the lower-dimensional vacua have 
positive vacuum energy, they are metastable, with local regions eventually transitioning back to the D-dimensional 
de Sitter space. This recycling process continues indefinitely, fragmenting future infinity into regions with diverse 
properties. 

When a compactification or interpolating solution is nucleated out of the D-dimensional de Sitter space, it occupies 
nearly a Hubble volume of the background space. This can be seen by noting that both the size of the q-sphere 
(through the location of the potential maximum Eq. [25]) and the amplitude of a (through uj~^ Eq. [76)1 is set by A. 
Thus, the naive picture that emerges is that a nucleation event corresponds to taking roughly one horizon volume of 
the background dS space and replacing it with the post-nucleation spacetime. In the asymptotic future, this patch 
then evolves into two regions: one in which the spacetime approaches a p + 2-dimensional vacuum, and the other in 
which it returns to the D-dimensional de Sitter space. 

At the very most, the region containing the p + 2-dimensional vacuum can correspond to removing a comoving 
Hubble volume from the background dS space. The interpolating solution is anisotropic when viewed from the full 
ZJ-dimensional theory, which introduces additional complications. Leaving this subtlety aside, we can course-grain, 
and attempt to determine how the comoving volume of the original de Sitter space is distributed between the various 
vacua. 

Consider a fixed parcel of comoving volume that is on some fixed initial value surface everywhere in the D- 
dimensional de Sitter space. Extending comoving geodesies from this initial value surface, we can trace the evolution 
as a function of proper time of the fraction of such worldlines that remain in the original vacuum Vdsi^) those 
that encounter a nucleation event 'Pi{T) of type i. If we restrict our attention to the volume cut out by the nucleation 
events, and do not attempt to follow geodesies into the interpolating region, this should be a valid procedure. This 
type of volume weighting is identical to the procedures applied in various measures for eternal inflation (see for 



example 31|, whose conventions we use below, for further discussion). 



Conservation of probability implies that 

Pd5(T)+^P.(r) = l. 

i 

Disregarding decompactification transitions back to the D-dimensional dS space, the rate equations are 

dVds 



dr 



Y.<''P'^s{rl (89) 



5 = >^fVds{r\ (90) 
where the Ki are the probability per unit proper time to experience a transition, given by 

Kf = Vol (r!z,_i) i7^-j^-^'r,. (91) 



The rates per unit D- volume Ti are given for the various types of transitions by Eq. 1801 and Hds is the Z3-dimensional 



30 



Hubble constant. On the initial value surface Vds{T = 0) = 1, and solving the rate equations we obtain 



Vds = 



^^dS 



(92) 
(93) 



The comoving volume fraction remaining in the Z?-dimensional dS space goes to zero (however, the physical volume 
will be exponentially increasing). At late times, the distribution of worldlines reaches an equilibrium state, with the 
relative fraction distributed between any two vacua i and j given by 



— — ^ exp 



^dS 



{a, - aj) 



(94) 



Therefore, under this measure, vacua with the largest possible values of a are exponentially favored. 

Taking into account the fact that there are decompactification transitions back to the D-dimensional dS space 
requires us to follow worldlines between the vacua of different dimensionality. This will undoubtedly involve a number 
of important subtleties, but we might hope that a naive addition to the rate equations provides a rough picture of 
what happens: 
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where the rates are given by Eq. [SHI The full solution to the rate equations can be found using the methods 
outlined in [33 |. but considering a single transition type, the solution is 
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Again, at late times, a steady state behavior is reached, with the relative weight of the two vacua given by 
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(99) 

(100) 



which is typically much greater than one, since jS'^s^l 3> jS'Jgl when the scale setting the p + 2 dimensional vacuum 
energy is much smaller than the scale setting the D-dimensional vacuum energy. This result is in accord with the 
qualitative expectation that a statistical system spends most of its time in the state of highest entropy, and was 
commented on originally in Refs. [1, Here, it comes about because the rate to evolve back from the p+2 

dimensional vacuum to Z)-dimensions becomes arbitrarily small as the p + 2-dimensional vacuum energy goes to zero. 

This is only one of many possible measures over the vacua produced from the background dS space. Connecting 
such measures to observable quantities is inherently ambiguous. For example, it is not obvious how to compare the 
number of observers in each unit of comoving volume, both because it is unclear what a unit of comoving volume 
has to do with observers, and because we are comparing infinities and there is no unique regulator (notwithstanding 
the fact that one must compare volumes in vacua of different effective dimensionality) . We have little to add to this 
discussion, but it should be possible to generalize many of the existing measure proposals to this scenario. Another 
possibility is offered by the no-boundary proposal. Here, each Hubble volume of the background dS space is treated 
as an independent system, "nucleated" with probability set by the various instantons found above. The relative 
probability of obtaining various solutions is proportional to the ratio of rates, as in Eq. 1941 

We also note that there might exist a progressive relaxation of the number of dimensions if the appropriate q-iorui 
fluxes were available. In this scenario, after compactifying to Mp+q2+2 x S"^^ j where Mp_|_q2+2 is a p-l- 92 + 2-dimensional 
vacuum solution, a second transition could produce a product compactification Mp+2 x 5"^^ x S'^^ . Including electric 
flux as well, one could consider the nucleation of charged membranes in the p + 2-dimensional vacua, embedding the 
standard four-dimensional picture of false- vacuum eternal inflation inside of the dynamical compactification scenario. 
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VII. THE COSMOLOGICAL CONSTANT PROBLEM 



Within the landscape of radion potentials, it is possible to find effectively 4-dimensional vacua with positive vacuum 
energy. In order to promote this toy model to a plausible theory of our universe, it is necessary to ensure that 
there is at least one vacuum with sufficiently small positive energy to account for the observed vacuum energy 
Pvac ^ 1-5 X 10^-^^^ [34]. This four-dimensional vacuum energy is determined by a cancellation between the curvature, 
flux, and higher-dimensional cosmological constant contributions to the radion potential Eq. 1201 There is one vacuum 
per value of the charge Q, and since the flux should obey a quantization condition Q = en with e the gauge coupling 
and n integer, there is a discrete number of vacua. 

The question of how miraculous the required cancelation is depends on how finely the vacuum energy scans with a 
change in the discrete values of the charge. If the gap between successive values of the vacuum energy is large, then 
the gauge coupling and/or higher dimensional cosmological constant would have to be extremely fine-tuned in order 
to produce the observed value of Pvac- However, if the gap is quite small, and there are many vacua with energies in 
the correct ball-park, then it is reasonable to believe that the correct vacua can be populated, even in the absence of 
fine-tuning. 

We find two ingredients in our scenario that reduce the amount of tuning in A and e necessary to naturally produce 
the correct value of the four-dimensional cosmological constant: a large number of extra dimensions q, and the presence 
of multiple quantized g-forms with incommensurate fundamental charges. In essence, these are the same methods 
used by Bousso and Polchinski [s^ to obtain vacua with a naturally small value of the cosmological constant, although 
here they arise in a different context. 

If there are multiple copies of a 5- form, then in the radion potential Eq. [201 is replaced by the sum 

Q'-t.Q"=i:e^ni (101) 

i=l i=l 

where is the gau ge coupling for each copy, and rii is the number of units of fundamental charge. Following Bousso 
and Polchinski [35j , we can visualize charge space as a J-dimensional grid of points with a spacing in each direction 
set by Bi. Inside of this charge space, construct a shell of radius Qo, and width AQ. In the bulk of this shell, the 
total charge produces a vacuum energy between approximately zero, which defines the radius Qq, and the observed 
value of the vacuum energy, which defines the width AQ. In order for it to be likely that there is a vacuum with 
the appropriate value of the cosmological constant, we require that the volume contained inside of the shell is greater 
than the volume of a unit cell in charge space 



'-j^ ^ Vol (a/- 1)^,7-1 



n 



< '—^Q-^-'AQ, (102) 



where d is a degeneracy in values of Vq. If all the Ci are different, then d = 2'^ due to the invariance under rii — > —rii 
(the total volume in the shell is reduced to the volume contained in a region with positive rii). It is possible that d 
can be somewhat larger if the are commensurate. 

Applying this analysis to the radion potential, the total amount of charge Qo necessary to produce a zero p + 2- 
dimensional vacuum energy can be found by substituting the approximate location of the minimum Eq. [21] into the 
radion potential, setting the latter to zero, and solving for Q. The result is an order of magnitude estimate for Qo 
given by 



[p + 1)9 V 2A 

For a given p + 2-dimensional vacuum energy pvac, we can define the width AQ by 

V dQ 



^Q = Pvac(^^^\Q=Qj , (104) 



where the term in parentheses is the derivative of the height of the potential minimum with respect to Q evaluated 
at Qo. Substituting the approximate position of the potential minimum given by Eq. [Hjinto the radion potential and 
evaluating the derivative, we can find the scaling with A 

-M';Uci,.,)(±.)''^. (105, 
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FIG. 21: A log plot of the function C{p, q) for p — 2 and various q, obtained numerically by finding the value of the potential 
minimum for values of Q that yield Vb ~ 0, constructing the function V{(f>-[Q]), evaluating the derivative, and factoring out 
the scaling with A determined in Eq. 11051 It can be seen that the behavior is approximately exponential in q. 



where C {p, q) is a function of geometrical factors only, and can be determined numerically. An example for p = 2 as 
a function of q is shown in Fig. 1211 where it can be seen that this factor decreases exponentially with increasing q. 
Re-arranging Eq. 11021 we require that 
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It is interesting and potentially significant that the effects of fine-tuning in D-dimensions are amplified in the lower- 
dimensional theory in a number of different ways. In particular, the higher-dimensional cosmological constant is 
raised to a power dependent upon the number of compact dimensions and the number of fluxes, the individual gauge 
coupling appear as a product, and there are potentially large geometrical factors. With many compact dimensions 
and/or fluxes, mildly tuned values A < M|, and < 1 are raised to large powers, vastly increasing the volume of the 
shell in charge space. 

For example, taking {e = .1,A = .1M|,} and q = 7, the inequality is satisfied for J > 18. With more tuning in A, 
it is possible to satisfy the inequality for smaller values of J and q. With A — lO^^M^, the inequality is satisfied for 
{q = 7,J > 5}, {q = 6, J > 6}, {g = 5, J > 8}, {g = 4, J > 10}, {q = 3,J > 15}, or {g = 2, J > 30}. The inequality 
can be satisfied for large numbers of extra dimensions or fluxes individually. For example, with q 60, one requires 
A ^ -^M^ and e ~ .1, while for q = 2, one requires J > 137 for the same e and A. Even for the minimal values 
{q = 2, J — 1}, after setting e ~ 1, the necessary tuning of A = 10~'**M|, is some 75 orders of magnitude better than 
if one had started with a four-dimensional theory. 

Although we will not explore the issue in detail, how finely the can lower-dimensional vacuum energy scans will 
also affect any potential statistical predictions of the cosmological constant. A key element in the original arguments 
of Weinberg [36| for an anthropic prediction of the cosmological constant is that there are many vacua inside of 
an "anthropic window" of acceptable vacuum energies. The validity of this assumption requires that the vacuum 
energy scans quite finely, and must be tested on a case-by-case basis. For example, an analysis of this issue in the 
context of the string theory landscape can be found in Ref. [13] . For a modern review of statistical predictions for the 
cosmological constant see e.g. (ssl, Isoj. 



VIII. INCLUDING AN EPOCH OF LOWER DIMENSIONAL INFLATION 



An analysis of the p + 2-dimensional cosmological constant relies on the properties of the p + 2-dimensional vacua 
alone. However, perhaps the most interesting feature of the solutions discussed in previous sections is the existence 
of a region containing a p + 2-dimensional FRW universe. The radion field dynamically evolves towards its vacuum, 
approaching at late times a p + 2-dimensional de Sitter, Minkowski, or big-crunch spacetime. In our universe, scalar 
field dynamics is thought to drive an epoch of slow-roll inflation, and it is a natural question to ask if we can embed 
this into our model. 

We will consider two separate models of inflation in this section, both of which require the addition of a scalar 
field Tp into the _D-dimensional theory. If this field has a positive mass squared, then it will be possible to construct 
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FIG. 22: Two models of inflation arising from the potential Eq. IllOl In both cases, the t/j field becomes destabilized between 
the minimum and maximum of the potential along the 4> direction (the profile along -0 = is shown as the solid line). This 
can either drive (left) or end (right) an epoch of slow-roll inflation. For the potential on the left, there exists an interpolating 
solution, which naturally sets the initial conditions near '0 = 0. Inflation then occurs as the fleld rolls in the direction. For 
the potential on the right, there is no interpolating solution, but the Euclidean compactiflcation solution can set the initial 
conditions for the field very near the potential maximum. Inflation occurs as the field rolls in the (j) direction, eventually ending 
when the field falls into one of the minima in the ip direction. 



the compactiflcation, interpolating, and Nariai solutions of Sec. |V] exactly as before. This is not a cosmologically 
interesting theory, but by coupling -0 to the D-dimensional Ricci scalar and field strength, it is possible to drive the 
mass squared negative for large enough curvature and/or flux. In the dimensionally reduced theory, this corresponds 
to couplings between the radion field (j) and the scalar with the large curvature regime corresponding to small (p. 

In the model shown in the left panel of Fig. [521 the instability causes if) to slowly-roU, while is stabilized in a 
minimum of the potential. In the model shown in the right panel of Fig 1221 (j) is the slowly-rolling inflaton, while if) 
is in an unstable equilibrium, eventually rolling off to a global minimum of the potential, and ending inflation. The 
slow-roll parameters in each case are satisfied near critical points of the potential, with slow-roll inflation occurring 
mostly in the ip direction for the potential on the left and mostly in the (p direction for the potential on the right. 

The coupling of '0 to curvature and flux also allows for a natural explanation of the initial conditions for inflation. 
For example, the potential on the left of Fig. [221 admits an interpolating solution. It is possible to adjust the couplings 
such that the instability in the ip direction develops only in the spacetime region containing the p -I- 2-dimensional 
FRW universe. In such models, slow-roll inflation in the -0 direction is triggered by the compactiflcation of the 
extra dimensions. The potential on the right of Fig. [22l does not admit an interpolating solution, but there is a 
compactification solution at the potential maximum along the <j) direction. Evolution away from the maximum in the 
(j) direction can drive slow-roll inflation, until the instability in the direction takes over, at which point inflation 
ends. 

We now examine the two models of inflation in detail. In both cases, the D-dimensional action we will consider is 
given by 



S = 



- 2A - 



2q\ « 



(107) 

The functions /(ip) and h{ip) specify the coupling to gravity and flux respectively, and ■0 has potential v{'ip). The 
properties of the ip field depend on background field and metric configurations through the couplings f{ip) and h{ip). 
It will be convenient to go to the Einstein frame in order to utilize the solutions found previously in Sec. |Vl After 
some manipulation, the action becomes 
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Performing the dimensional reduction as before and going to the Einstein frame of the 4-dimensional theory we obtain 
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the action 



where 
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For the moment, we will find it convenient to work with the field R instead of the canonically normalized radion field 



Both models are specified by: 
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where m has mass dimension q + 2 and A has mass dimension q + 4 — n. Assuming an exact shift symmetry ?/' - 
this corresponds to the case of minimal coupling of the ^ field to gravity and flux. 

The key feature we will exploit is the i?-dependent mass of the ip field. Examining the second derivative of the 
potential with respect to V' at "0 = and fixed R, we obtain 
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a^<0, then 



where m = m/M'j^'^ . The ip field is stabilized at large R as long as rri > £,i^l^jp-- However, if ^2 s,i;^+2 
-0 becomes destabilized at small R. The region where ijj becomes tachyonic allows for the possibility of an altered 
vacuum structure where tp can obtain a VEV. We assume that the second term in the potential V{ip) does not affect 
the properties of the potential in the vicinity of -0 — 0, but it will be important in establishing the vacuum structure. 

In order for the solutions described in Sec. |V]for the radion field to exist, we require in both models that the sign 
change in the second derivative Eq. lll2l occurs somewhere between the maximum and minimum of the potential along 
^ = 0: 
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(Af|,/A)(9-i)/2 and 



^ must therefore be of order kjMl for the sign change in the second derivative to occur in 



From the discussion in Sec. Ill Bl if the minimum along = is positive, then Q 
R-j- ~ A~^/2 The ratio — rr,-,-,L?+ +v.Qr.i3f,-ir-Q r^^Aar A / A r2 

the desired position. 

Consider first potentials of the type shown in the left panel of Fig. [22l where plays the role of the inflaton. For 
potentials of this type, the interpolating solution exists when the condition on the location of the second derivative 
Eg. IllBl is satisfied. The field configuration on the big-bang surface of the 4-dimensional open FRW universe inside of 
the interpolating solution is set by the value of ~ 0_ at the stationary point of the interpolating solution and the 
stabilized value ~ 0. Near the minimum of 0, Eg. II 121 becomes negative and ip is destabilized and begins to slowly 
roU. 

These are essentially "hilltop" models [40|]; near the minimum in 0, the potential is very close to quadratic over 
a range of ip of order M4, after which the potential steepens, and inflation ends. As is near its minimum, we can 
analyze this as a single-field model of inflation, where tp is the slowly rolling fleld. Inflation occurs near the maximum 
of the potential in the "0 direction which this guarantees that the flrst slow-roll parameter 
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is small, but not the second slow-roll parameter rj 



(115) 
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FIG. 23: Numerical evolution in the presence of the potential in the left panel of Fig. 1221 The parameters specifying the 
potential are {g = 2,n = 6, Q = 142, A = 5 x 10"^m = 3.24 x 10"'', = 1.63 x 10"^, ,^2 = -9 x 10"^A = 3.3 x 10"^^}. The 
stationary point of the interpolating solution sets the initial condition for i/) to <^o = 9.6. We choose t/; = .3 to obtain enough 
inflation (for this initial condition, there are approximately 110 e-folds). 



When both are small, we will find models with greater than 60 e-folds of inflation. If e is negligible, the spectral index 
is 

n, -l~-2|r?|. (116) 

Using the WMAP 5-year [s^l central value of Ug = .96, we require \ri\ = .02. We must also ensure that the scalar 
power 

is of the correct magnitude Vn ~ 2.5 x 10~^. This, together with the absence of tensor perturbations in the CMB, 
constrains the scale of inflation to be V^^^ < 2 x 10^^ GeV. 

Satisfying the constraint on the spectral index requires m, ^i, an ^2 to be somewhat less than one so that the 
second derivative Eq. 11121 is sufficiently small. This represents a true tuning, and corrections will introduce the so- 
called "77-problem" (corrections drive 77 ~ 1). The correct scalar power can be obtained by a combination of initial 
conditions and either taking A/Af|, ^ 1 (scaling the entire potential) or choosing an appropriate value of Q such that 
the minimum along ■0 = is of the appropriate height. One model that satisfies all of the observational constraints 
has the parameters {Q = 142, A = 5 x 10"^ m = 3.24 x 10""^, 6 = 1-63 x 10"^, ^2 = -9 x 10"^}. An additional 
tuning of A is necessary in order to end inflation in a vacuum with very small cosmological constant. In this example, 
it is necessary to take n = 6 and A = 3.3 x 10""'^^. 

The numerical evolution of the fields in this model is shown in Fig. 1231 where we have transformed from R to the 
canonical radion field 0. The initial condition for the field is set by the stationary point of the interpolating solution 
which is approximately at ~ 9.6. Because the ip field is stabilized in the _D-dimensional de Sitter space, it is deposited 
at -0 ~ modulo fluctuations. We might expect the size of typical fluctuations to be of order Aip H ^ 10"^ in this 
example. The universe begins dominated by curvature, and because the gradient in the ip direction is quite small, 
the field moves in the direction, and begins damped oscillations about the minimum in cj). This begins an epoch of 
inflation, as the field slowly rolls in the ip direction, as can be seen in Fig. [221 The total number of e-folds depends 
on the initial condition for ip, which in this example is taken to be ipQ — .3. This is much larger than the expected 
size of fluctuations, and so the total number of e-folds will typically be vastly larger than the required 60. In fact, 
the potential in the ip direction is sufficiently flat that an epoch of slow-roll eternal inflation can be triggered inside 
of the interpolating solution. As long as there are enough e-folds (around 60, which is not surprisingly the number 
typically needed to solve the flatness problem), the original curvature of the FRW region in the interpolating solution 
is diluted below the current bound from observations of the CMB. This is identical to the dilution required in the 
epoch of "open inflation" that occurs inside of bubble universes, see e.g. [41l] . 

We now move to the second type of inflationary potential, shown in the right panel of Fig. 1221 where (j) plays the 
role of the inflaton. There is no interpolating solution for these models, and the initial conditions are set by the 
compactification solution. This places the field in the near-vicinity of the potential maximum in the (j) and close to 
the stabilized value i/j = 0. If the field starts close enough to the maximum, an epoch of slow-roll eternal inflation 
will occur. Spacetime regions fluctuate away from the maximum, beginning slow-roll inflation with a variety of initial 
conditions. 
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FIG. 24: Numerical evolution for the potential in the right panel of Fig. 1221 The parameters specifying the potential are 
{g = 2, n = 6, Q = 3651.5, A = 1 x 10"'^, m = 1.32 x 10"'^, = .199, ,^2 = -1, A = 3.6 x lO"'*}. The Euclidean compactification 
solution places the field near the maximum of the potential in the 0-direction near tp = 0. We choose the initial conditions to 
be {(f>o = 16.122, -(/jo ~ .003}, which yields approximately 90 e-folds of inflation. 

This model is similar to hybrid inflation [13]; as the field slowly rolls in the cj) direction, an instability develops in 
the tjj direction that ends inflation. In these models Q ~ Qmax, and an approximate inflection point develops in the 
(j> direction. The scale of inflation is therefore set exclusively by the overaU height of the potential, and it is necessary 
to take A/M|, ^ 1 to satisfy the constraint on the scalar power. However, because we no longer need to tune the 
second derivative in the "0 direction to be small, we need only ensure that V' becomes destabilized at an appropriate 
value of (f). This requires to ^ A/Af|, with and 1^2 order unity. For example, a working set of parameters is 
{Q = 3651.5, A = 1 X 10"^ to = 1.32 x 10"^ = .199, 6 = -!}• Evolution from {^o = 16.122, t/j = .003} yields the 
trajectory shown in Fig. [211 and produces a spectrum of fluctuations consistent with observation. Again, an additional 
tuning of A is necessary to end inflation in a minimum with small vacuum energy. In this example, we must take 
n = 10 and A = 3.6 x 10"^. 



IX. CONCLUSIONS 

Despite the enormous amount of research into theories of extra dimensions, very little is known about the possible 
ways to make a transition from a higher-dimensional to lower-dimensional universe. In this work we have explored 
how such a transition might occur in one particular setting, and furthermore explored the landscape of possible vacua 
of different dimensions and vacuum energy that would emerge. 

The starting point for this dynamical compactification scenario was D = p + q + 2-dimensional Einstein gravity in 
the presence of a positive cosmological constant and a set of g-form field strengths. In this theory, the Z)-dimensional 
de Sitter vacuum solution is semi-classically unstable to the nucleation non-singular interpolating and compactification 
solutions that contain a region in which q of the dimensions are compact and there is a p -I- 2-dimensional cosmological 
spacetime. This region evolves towards either a big-crunch singularity or a p -I- 2-dimensional universe with zero or 
positive vacuum energy, perhaps undergoing an epoch of slow-roll inflation in the process. In this way, it is plausible 
that our universe arose from such a dynamical compactification of some number of extra dimensions. 

The properties of the p + 2-dimensional vacua depend not only on fixed quantities such as the higher dimensional 
cosmological constant and total dimensionality, but also on variable quantities such as the charge and properties of 
other matter fields. Since the charge can take a number of values, and there can be multiple kinds of flux, a landscape 
of vacua with different vacuum energy and effective dimensionality emerges. Each of these vacua is populated by 
the dynamical compactification process described above. This mechanism of landscape population is fundamentally 
different from false vacuum eternal infiation, which operates entirely within the realm of the lower-dimensional effective 
theory. Since the interpolating solutions are close cousins of non-extremal p-branes, their semi-classical production 
is more analogous to the nucleation of charged black holes in dc Sitter space. 

In any theory where a number of vacua are realized in different spatiotemporal regions, the prediction of the 
parameters in a low-energy effective theory becomes statistical at best. Such a prediction will necessarily involve 
both dynamics and selection effects. Before tackling such complicated questions, a zeroth order requirement is for 



In the landscape of Type IIB string theory, the volume modulus can provide an obstacle to populating vacua through eternal inflation [43|| , 
and it is interesting to identify an alternative mechanism. 
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the theory to contain at least one vacuum with the observed properties of our universe. In our toy model, we can 
ask if there are vacua with the correct dimensionality and the observed value of the cosmological constant. Finding 
vacua of the correct dimensionality requires that there exists a {D — 4)-form field strength. Naturally obtaining the 
observed value of the cosmological constant requires that there is 1) tuning of the higher dimensional cosmological 
constant, 2) tuning of the fundamental unit of charge associated with the flux, 3) multiple types of g-form flux with 
incommensurate charges, 4) many extra dimensions, or some combination thereof. This is very similar to the story 
in the string theory landscape. 

To gain some understanding of the dynamics, we have explored the phenomenology of nucleation rates. Generically, 
the rates are suppressed by some appreciable fraction of the higher dimensional de Sitter action, and are higher for 
smaller values of the lower dimensional vacuum energy (the highest rate is to the most negative minimum). Another 
interesting behavior arises when comparing the rate to vacua of different dimensionality, which at fixed normalized 
charge Q/Qmax is maximized for a particular p or range of p depending on the total dimensionality. Comparing the 
rate into a p + 2-dimensional vacuum to the rate for decompactification back to D-dimensions, we find in agreement 
with previous studies that Minkowski vacua are completely stable and the rate into a p + 2-dimensional vacuum with 
a small positive vacuum energy is much higher than the rate out. These rates can be used to examine the structure 
of the higher-dimensional multiverse, for example calculating the distribution of comoving volume. 

Although we have presented a relatively complete picture of the solutions and their nucleation, we have neglected 
a number of potentially important effects. Perhaps the most pressing is the question of stability. It is well known 
that black branes in flat space are unstable [44], which may have implications for the stability of the interpolating 
solutions. The stability of the dSp+2 x S'^ compactification solutions was examined in Ref. (sl] (see also [H, |46|). 
where an instability was discovered for g > 4. This is perhaps troubling since the existence of more than 3 extra 
dimensions is theoretically well motivated from string theory. 

Recently this problem was revisited in Ref. [47j . who found stable warped solutions for 9 > 4, and proposed that 
such solutions are the endpoint of evolution from the unstable solutions without warping. In this case, the existence 
of an instability would not affect the qualitative aspects of our analysis in any number of dimensions, but addressing 
this question in more detail is clearly necessary. 

In addition, there are many are many directions in which one could proceed to make more detailed and quantitative 
predictions in this multiverse model. This necessarily involves determining an appropriate measure over the various 
solutions, but there may be other more direct signatures. For example, in false vacuum eternal inflation, it is known 
that the inevitable collisions between bubble universes might lead to direct signatures of other bubble universes [i^ . In 
complete analogy, we might expect that some of the interpolating solutions overlap, leading to potentially observable 
inhomogeneities in the p + 2-dimensional cosmologies. This situation is also analogous to multi-centered black hole 
and black brane solutions. 

It is hoped that this model will provide different insights into issues surrounding extra dimensions and theories with 
multiple vacua than those obtained by studying eternal inflation in the string theory landscape. 
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A more general p-l- 2-dimensional metric ansatz is to assume homogeneity but not isotropy. In four dimensions, these 
are the Bianchi-type cosmologies. Here, we detail the construction of solutions without spatial curvature, although 
it is possible to generalize our methods to all homogenous solutions in any number of dimensions. For flat spatial 
sections, the most general homogenous metric can be written as 
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APPENDIX A: HOMOGENOUS BUT ANISOTROPIC SOLUTIONS 



P+i 




(Al) 



where the Xi are cartesian coordinates that generally range between — 00 < a;^ < 00. There are p + 2 (dimensionless) 
scale factors ai(r) which are all arbitrary functions of t. 
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It is possible to decouple the Einstein equations by finding equations of motion for the field and the product of 
metric functions 

P+i 

C=Y[a,. (A2) 



i=l 



Making this substitution yields 




C_2(p+1) V 



C P Ml 



(A4) 



p+2 



''i^+-^4>^-Ml-^^V'. (A5) 

The constant E parametrizes the degree of anisotropy, and is related to the metric functions ai (r) as we outline below. 

The Ricci scalar is again given by Eq. [35] (as can be explicitly checked using the equations of motion) , and so as long 
as the field and its derivatives remain finite, the C = surface is completely non-singular. From the field equation 
Eq. IA5[ this requires that 0(r = 0) = 0. The C, = surfaces correspond to event horizons in the D-dimensional 
geometry, just as they did in the black hole example of the Introduction. From Eq. IA3[ if the field and its derivatives 
stay finite, then the anisotropy term always dominates at small C,. This leads to the universal behavior C, — Sr as 
T ^ 0, which can be recognized as the Kasner solution. 

The value of S and the time dependence of the individual Oi(r) can be found from the remaining Einstein equations. 
Near r = 0, the metric is always in the form of the Kasner solution, where the scale factors are a power law in r 

a,{T)^af\^^. (A6) 

There are two Kasner conditions 



P+i p+i 



i=l i=l 



that constrain the exponents ni. It can be shown that unless there is a single non-zero n^, at least one of the must 
be negative (see e.g. Ref. E^). This implies that if there is more than one non-zero rii, at least one scale factor 
diverges at r = 0, while the others shrink to zero. If it is possible to continue across the t — Q surface to a region 
where r is spacelike, only one of the metric functions can go to zero (otherwise the analytic continuation becomes 
ambiguous). Therefore, in the vicinity of r = 0, only one of the Kasner exponents can be positive. From the Kasner 
conditions, in order to obtain C = with a single positive Kasner exponent, all of the other exponents must be zero. 
From the equations of motion, the anisotropy parameter E is given by 

S^lfaf). (AS) 

i=l 

The full time dependence of the scale factors can be determined by identifying a number of constants of the motion. 
For example, specializing to p -I- 2 = 4, and following Ref. [s^l, there are three constants of the motion 

111 = ai (0203 - 0302) , 

112 = 02 (0301 - 0103) , 

113 = 03 (0102 - 0201) . (A9) 



One then obtains equations for ai(T) and 02 (r) 



di 


IC 


1 


ai 


3C 


3C 


0.2 




1 


0.2 




3C 



(n2 + 2ni) , (Aio) 

(Hi + 2U2) . (All) 
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The anisotropy parameter E is given in terms of the constants of motion by 

Y,'^ = Uj+Ill+UiU2. (A12) 

The time-dependence of 03 can then be determined from (. This procedure can be appUed in an arbitrary number of 
dimensions by again identifying the appropriate constants of the motion. 



1. A = 



We now construct anisotropic solutions in the absence of a D-dimensional cosmological constant. Before analyzing 
the dimensionally reduced theory, we briefly discuss some analytic solutions that can be generated from an anisotropic 
metric ansatz. Included in this set of solutions are the extremal and non-extremal p-brane solutions of Ref. [51,] . The 
metric for such solutions in D-dimensions is given by [l5| 



1 - 



R+ 
B 

R+ 
R 



9-1 



1 - 



9-1 



-, 1^ 



1 - 



~R 



q-l 



dx ■ dx 



R- 



dR^ 



where i?+ and i?_ are the two event horizons related to the charge by 

p2 _ (g+P)(g- 1) ,r> R 

^ - 2{p+l) (^+^-) 



(A13) 



(A14) 



There is a curvature singularity at i? = i?_ . When R^ = R_ , the event horizons become degenerate, and the solution 
is said to be extremal. These are none other than the extremal solutions we constructed in Sec. IIVI It is also possible 
to have super-extremal solutions, in which there is a naked singularity (these are similar to the timelike singular 
solutions of Sec. IIV|) . 

In the region between the horizons i?_ < R < i? is a timelike variable. Near the horizon at we now 
show that the metric is of the Kasner form. This can be done analytically ioi p — 2 by dimensionally reducing to 
4-dimensions and performing a conformal transformation to the 4-dimensional Einstein frame. The metric Eq. IA13I 
becomes 



ds-" 



{MoRf 



+ {MdRY 
Changing coordinates to 



R+ 
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R+ 
R 



9-1 



R- 
~R 



dR^ + {MdR'') 
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- 1 
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~R 
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dt'. 



L — Rj^ — R, 

and expanding to lowest order around the horizon at L = 



[MdR+Y i?4 



ds^ = - 



9-1 



i?4 



(9-1) 

(g - 1) {MpR+Y 

q-l^ 1/3 
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(MdR^ 



Ldr 



R+ 



Re-scaling L, x, and t by the constant factors, we obtain 



ds^^-±dL^ 
L 



1 - 



dx ■ dx ^ dt^ . 

■r+ 



R- 
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q-l 



dx ■ dx 
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(A15) 
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Finally, defining the proper time variable r = we recover the metric 

ds^ = -dr'^ +dx-dx + -^df. (A19) 
r+ 

This is of the Kasner form, with all but one of the exponents zero as expected. 

The non-extremal p-brane solutions discussed above can also be obtained from the dimensionally reduced theory by 
evolving in the presence of the radion potential, just as in the analysis of Sections IIIIl IIVI andfVl The non-extremal, 
extremal, and super-extremal p-brane solutions correspond to different initial conditions for the radion field. In fact, 
these are the only solutions that possess an asymptotically fiat D-dimensional region, and the radion and scale factor 
profiles will be qualitatively similar to the extremal and singular (timelike and spacelike) solutions found in Sec. IIVI 

Our method of constructing solutions in the dimensionally reduced theory will be identical to the procedure outlined 
in Sec. IIIIl where spacetime regions with timelike and spacelike r are matched across event horizons. We will discuss 
the evolution of the radion field in each case separately. 

Returning to the equations of motion for C, it can be seen from Eg. IA4I that the sign of the second time derivative 
of C is fixed by the sign of V . For A = 0, the potential for is negative in the vicinity of the potential minimum 
and approaches zero from below as (/>—!■ oo. Evolution to — > oo is therefore in a spacetime region where r is a 
spacelike variable and V —V. Just as for the maximally symmetric case, there are attractor solutions approached 
from generic initial conditions. The analysis is similar to that of Sec. IIIIl 

We now look for solutions in regions where r is a timelike variable as in Eq. lAll Concentrating on motion in the 
vicinity of the potential minimum, we can apply the methods used in Sec. IIIIl Since the minimum is negative, the 
scale factor is bounded, with the frequency of motion given from Eq. IA4I bv 

co',^^^^\V{^-)\. (A20) 

Substituting into the field equation Eg. IA5I the solutions for (f> can be written in terms of Gegenbauer polynomials, 
with the index determined by 

which, substituting with Eq. 1581 yields a = 1 independent of p. Therefore, for solutions with stationary points 
infinitely close to the minimum, there are solutions that interpolate across the minimum exactly once, and the 
oscillatory solutions found in Sec. |V] cannot not exist. For an initial condition displaced slightly from the minimum, 
the frequency increases, causing the solutions to go singular in the vicinity of the would-be second stationary point. 

Putting the full solution together, if there is a non-singular stationary point on the right side of the potential 
minimum, then it is possible to connect the asymptotic region — > oo with r spacelike across an event horizon to 
a region in the basin of attraction of the minimum with r timelike that contains a singularity. This is exactly the 
expected behavior for the non-extremal p-branes discussed previously. The causal structure will resemble the singular 
(spacelike) solutions of Sec. |Vl If there is no non-singular stationary point, then the geometry will contain a timelike 
naked singularity. The causal structure in this case will be identical to the singular (timelike) solution described 
in Sec. |Vl The extremal solution will be identical to the solution obtained in Sec. IIVI since the regime where the 
anisotropy dominates the energy density is pushed to r ^ — oo. 

2. A > 

When A > 0, the potential approaches zero from above, implying that r is a timelike variable in the asymptotically 
il'-dimensional region of the geometry. There are attractor solutions where r is timelike interpolating between a 
non-singular stationary point and </> oo of the type discussed in Sec. [V] Continuing across the stationary point to 
a region where r is spacelike, we can look for motion in the vicinity of the inverted potential maximum. The solution 
for the scale factor is oscillatory and the solution for the field is again a Gegenbauer polynomial with index given by 
the solution to 

a^+a= / 1^^^)'. (A22) 
2(p+l) |y(0+)| ^ ' 

The solution for a is always less than one, implying that there are no non-singular trajectories crossing the maximum 
of the potential - the field cannot evolve quickly enough to traverse the maximum before the scale factor hits its 
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second zero. Thus, in contrast to the solutions in Sec. IIVI for an FRW metric ansatz, there can be no interpolating 
solutions connecting the asymptotic (j) —^ oo region to the potential minimum. The geometries will therefore have 
a causal structure similar to the singular solutions of Sec. IIVI Moving the stationary point in the motion of the 
radion field to the potential maximum, the anisotropy becomes unimportant, and it will be possible to find the Nariai 
solution discussed in Sec. [IV] 
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